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LONGITUDINAL STABILITY AND CONTROL 
OF A DUCTED-ROTOR FLYING PLATFORM 


INTRODUCTION 

The ducted-rotor Flying Platform is essentially a stand-on type, VTOL, 
lifting device operated by a pilot. Lift is produced by two, counter-rotating, 
non-articulated, rotors located below the pilot's stand-on platform and 
powered by reciprocating engines. (Fig. 1) 

The conventional control systems of airplanes and helicopters are 
eliminated, since control is primarily obtained by the method often refer¬ 
red to as "kinesthetic control". (Ref. 1) 

Longitudinal and lateral control is maintained when the pilot, acting 
as an automatic feedback by his instinctive balancing reactions, provides 
the correcting moment necessary to control and stabilize the machine after 
a disturbance. In steady forward flight he provides the required control 

moment by shifting his weight and leaning in the direction of the desired 

£ 

motion. 

Directional control is obtained through the movement of yaw vanes 
located in the duct and operated by the twist grip handle bars, which also 
provide a method for power adjustment. 

In this Report, all reference? to the term ' Flying Platform" snould 
be construed to mean the subject vehicle. 

The purpose of this Report is f o predict and understand, by 
theoretical and ..^iytical studies, the attainable flying qualities asso¬ 
ciated with the stability and control of the ducted-rotor Flying Platform 
as affected by changes in duct design, vertical canter of gravity location 
and artificial stabilization, 4 - 







2 

The ducted-rotor Flying Platform has been flying for many years 
and several models have been built. 

Stability and control problems have developed for the longitudinal 
mode of motion of these machines. The main difficulties are associated 
with large nose-up pitching moments, poor gust response, large tilt 
angles for forward flight and marginal handling qualities. 

For the studies in this Report, equations of motion based on small 
perturbation theory were developed and the static stability derivatives 
were calculated from data obtained by wind tunnel tests performed at 
the David Taylor Model Basin. (Ref. 2) 

Variations in vertical center of gravity location, together with radi 
cally different duct designs, (Fig. 2), were analyzed for their effect on 
the stability and control of the vehicle. 

The effects of artificial stabilization of the type used on vehicles 
of this design were also considered. 
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NOTATION 


a 

L 

c 

f 


k, 

1 


/ 01 ' 


rn 

P 


q 

t 



z 


z^ss* l p cos a 


A, B, C / D, E 

a h' B h' c h' d h 


slope of blade section lift curve per radian 
number of blades 

o 

chord length of duct (ft.) 
equivalent flat plate area of pilot 
(approximately 9 sq. ft.) 

non-dimensional moment of inertia parameter (0.0108) 
radius of gyration (ft.) 

distance between drag center of pilot and center 
of gravity of gross weight aircraft (ft.) 
mass of complete aircraft (slugs) * 
poles for root locus for C 

m 

P 

dynamic pressure (lbs. per sq. ft.) 
time (sec.) 

distance between the of action of thrust force and 
center of gravity of gross weight aircraft (ft.) 
zero for root locus for C 

m 

distance between line of action of pilot drag force 
and center of gravity of gross weight aircraft 
(ft.) 

distance between line of action of H-force and 
center of gravity of gross weight aircraft (ft .) 
coefficients of the characteristic equation in 
forward flight 

coefficients of the characteristic equation in 


hover 
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increments to coefficients of characteristic 

equations for C and C in root locus equation 
mm 

M * 

drag coefficient of pilot where C ' = f/S 

D R 

H-force coefficient 
pitching moment coefficient 



thrust coefficient 
weight coefficient 


static stability force and moment 
coefficient derivatives, (moments 
taken about the quarter chord) 


overall values of force and moment coefficient 
derivatives 


D = fp V / 2 drag force of the pilot (lbs.) 

horizontal force (lbs) 

F X total force measured along X body axis (lbs.) 
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F 


Z 


H 

Iy = mk y 2 


It 

L.C.D. 

M = M., - M 


M 


c/4 



R 

£lBP 

v- ! 

S.C.D.. 

T 


V 

w 


total force measured along Z body axis (lbs.) 

total force parallel to rotor disc (lbs.) 

pitching moment o t inertia of gross weight 

2 . 

aircraft (approximately 125 slug-ft. ) 
gain of artificial stabilization feedback loop 
Long Chord Duct 

pitching moment excluding pilot control moment, 
measured about Y axis (ft.-lbs.) 
aerodynamic moment created by the aircraft, 
without pilot, measured about the c/4 
point, (ft.-lbs.) 

control moment exerted by the pilot (ft.-lbs.) 
total moment measured about the Y axis (ft.-lbs.) 

period of oscillation (sec ) 

also poles for root locus for C 

a 

rotor radius (3.5 ft.) 

also Routh's discriminant 

abbreviation for Root Locus Break Points 
2 

area of rotor (ft. ) 

Short Chord Duct 

total thrust, force of duct and rotor measured 
normal to the rotor disc (lbs.) 

*lso time to half or double amplitude (sec.) 
also time constant for phase lag of artificial 
stabilization feedback loop 
resultant velocity of aircraft (ft. per sec.) 
gross weight of complete aircraft including 
pilot (550 lbs.) 




X,Y,Z 

Z 

a 


7 

e 

x 

V 

P 

m 

T= piTR^TR 
ft 


General 



A( ) 


6 

body a^is system 
zero for root locus for C 

m 

a 

angle of attack, angle between relative^ind 
vector and plane perpendicular to the rotor 
axis (radians) 

flight path angle measured with reference to 
horizon (radians) 

aircraft pitch angle with reference to horizon 
(radians) 

root of the characteristic equation 
tip speed ratio or velocity ratio 

air density (slugs per cubic ft.) 
aerodynamic time {pec.) 

rotor angular velocity (radians per sec.) 

» * 

differential operator forms for time ratio t/r 

differential operator forms for time t 

subscript zero designates initial condition 
before perturbation 
designates perturbation quantities 
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DEVELOPMENT 

Equations of Longitudinal Motion 

The development of the equations of motion for this system Is based 
on the commonly used methods for the study of aeronautical dynamics. 

The basic theory is the familiar small perturbation theory that assumes 
the disturbed motion of the aircraft to be one of small oscillations about 
some steady-state flight condition. 

As usual, the equations of motion are written in accordance with 
the Newtonian laws of motion with the important requirement that all 
accelerations and rates of change of moments of momentum must be 
expressed along alkes fixed in space. 

Equations of motion, as well as other important work in this 
area, have been developed at the Advanced Research Division of Hiller 
Aircraft Company by G. J. S^sslngh and A. H. Sacks. 

Also, dynamic stability of a machine similar to a Flying Platform 
has been treated at Princeton University in Refs. 3 and 4. 

The development of the equations of motion in this Report follows 
the latter work, the principle differences being in definition and notation. 

For the specific case of the longitudinal motion of the Flying 
Platform, we have a dynamic system with three degrees of freedom; namely, 
translation forward or backward, translation upward or downward and 
rotation about the pitch axis. 

This system will be described sufficiently by three equations of 
motion for the so-called control-fixed case similar to that studied in 
airplane dynamics. The control-fixed case of the Flying Platform can 
be thought of as the condition where the pilot remains rigid on the 
Platform, from the initial condition, and induces no new or additional 
increments of control movement throughout the motion to be studied. 
















8 

In order to facilitate the description of these motions, choice is 
made of the Eulerian or moving axis system fixed within the aircraft and 
always moving with it. For this axis system, the standard N.A.S.A. 
right-hand axis system is selected with the positive X, Y, Z, axes 
in the forwaro, right side, and downward directions respectively. 
Choosing these axes as body axes, the origin will always be located at 
the center of gravity and will move slightly with respect to the aircraft 
whenever the center of gravity is shifted. The positive X axis is directed 
toward the front of the aircraft, parallel to the rotor plane, and 
aligned with the relative wind in azimuth but not necessarily in pitch. 

The positive Z axis is directed downward, perpendicular to the rotor 
plane, i.e. aligned parallel to the thrust line or rotor shaft. 

The Y axis will be the aircraft pitching axis located at the 
center of gravity. This axis system, with notation for the forces, 
moments, and displacements, is shown in Fig. 3. 

In this discussion, the forces and moments to be resolved along 
and about the body axis system, together with other notation, are 
described as follows: 

total moment measured about the Y axis located at the 
center of gravity. Clockwise or nose-up moments 
are considered positive. 

aerodynamic moment created by the aircraft without pilot 
and measured positive clockwis&-~(nose-up) about 
the quarter chord point. (This notation arises from 
the wind tunnel balance design for measuring H- 
force and moments in connection with data re¬ 
duction when transferring moments to an arbitrary 







c.g. 


M 


c/4 







axis located at the so-called quarter chord point. 

This quarter chord axis is parallel to the Y axis 
and passes through the thrust line at the quarter 
chord point.) 

total gross weight of the aircraft, including pilot and 
all components „ 

force measured positive upward along rotor shaft 

line, perpendicular to the rotor plane and parallel, 
but directed opposite to the positive Z direction, 
force measured perpendicularly to T force and in the X 
direction. It is taken as acting in the quarter 
chord plane parallel to the rotor plane. The H- 
force is positive (as in the case of forward flight) 
when acting aft or in the negative X direction, 
drag of the pilot measured parallel to the relative wind 
and acting at the di ig center of the pilot. Drag 
Center and center of gravity location of the pilot 
alone were assumed to coincide. Also, aerodynamic 
lift caused by the pilot was assumed to be neg¬ 
ligible . 

moment arm for T force or displacement of center of gravity 
from thrust line measured along X direction. Dis¬ 
tance is positive when center of gravity is displaced 
aft of thrust line, and negative when forward of 
thrust line as in the case of forward flight, 
moment arm for drag force of pilot where 1^ cosa 

and 1 is the distance from the total center of 
gravity position to the drag center of the pilot 


as shown in Fig, 1 . 
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2 moment arm tor H-force or distance along Z 

ri 

direction, between the center of gravity and 
the line of H-force. 

a angle of attack, measured between the rotor tip path 

plane and the relative wind where nose-up is 
positive alpha. (Alpha is negative in forward flight.) 

0 attitude or pitch angle, measured between rotor Up path 

plane and the horizon where nose-up is positive 
tt eta. 

7 flight path angle measured to the horizon. As shown in 

Fin. 3, 7 = 0 - a where gamma is negative as shown. 

In the development Df the equations of motion of the Flying Platform, 
it is assumed that the disturbed motion of the aircraft is one of small 
displacements or oscillations about each steady-state equilibrium-trim 
flight condition. This steady-state equilibrium-trim condition is achieved 
whenever the summation of all the aerodynamic and gravity forces acting 
on the aircraft equals zero. Also, with regard to the inertia terms, the 
Newtonian laws of motion apply to any axis system, provided that the 
forces and moments are resolved along and about the particular axis 
system and that the accelerations expressed along and about these 
axes are the true accelerations with reference to fixed space. 

With these considerations the dynamic equations of motion may be 
written in the form: 


£F = E 


+ F + F 

'x_ 'x 


- 0 


/4£RODYNAM)C v GPAV/TY 


INERTIA 


EF = F 


+ h_ 


^AERCL YNAMIC ^GRA.v' T y 


+ F 


-o 


INEPT/A 


EM, =i1 


■h 


Yc.q. 


M 


o 


V 


AERCDYNAA//C 


c.g. 


INERT,M 


(la) 

(lb) 

(lc) 




il 


Resolving the aerodynamic and gravity contributions along body 
axes (Fig. 3), the steady-state equations of. motion become: 

E + E r ~ H ~ Dccso^- 6 = 0 (2a) 

X A X S 

|T *h F = -J- D Sin oC +\a/ COS 0 ~ 0 (2b) 

A G 

r^’^-H-.„*T-x T +D- s -o pc) 

Using the small angle assumptions where: 

SlTlof ~ c/" CGS &T « / 

(3) 

SL71 Q ~ @ Q “ZZ / 

the steady-state equations are: 

F x ^F x = - H - D - W @ - 0 (4a) 

p£ t F ~ 0 (4b) 

A 

M Vc.g^ = M c/+ -H'Z h +T'X t + D'* d "0 (4c) 

Expressing the inertia contributions by using the Newtonian 
frame of reference of fixed axes and body axes as the moving axes , 
the inertia terms become: 

# • * | 

F, - - /r/t, V c os of ~/rn V Y $i*i = - /rr V - 

(5a) 

IF. r 1/ Y css of-sw Vsin d ~ —/W ( 5b ) 
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My = — I 0 

C *1 / (50) 

These inertia terms arise when an aircraft accelerates, along a 

curved flight path as shown in Fig. 4, where the forces and moments are 
expressed as d'Alembert's inertia-resistance reactions. It should be noted 
that the 2 component of m$ namely, the contribution - m$ a, is small 9 

compared to the contribution of the Z component of centrifugal force mV? and 
is hereafter omitted from F g , since it is assumed to be a higher order term. 

No similar conclusions can be made about the relative size of the X com¬ 
ponents of these forces so that both contributions must be retained in F. 

Summing all the forces and moments as stated in Eq. 1, the 
longitudinal dynamic equations of motion for the Flying Platform become: 
H-Force Equation 


V 


(5a) 


(6b) 


EFx = -H-D-W0 -s** VY — o 

T iifust Equation 

£F Z = ~T-c<fD + W+ 'tjiVy = o 

Moment Equation 

£ M y,., = fV H ■ z H + T* x T +D • z D - 1 Y e - o (6c) 

In addition to the small perturbation assumptions, it is assumed that 
the changes in the external forces and moments acting on the Platform, 
because of the small departure from steady-state motion, will depend 
entirely on the displacement and disturbance velocities along and about 
the aircraft axis, and also that the frequency of the motions is low so that 
the external forces and moments do not depend on the accelerations involved 



*- ■»» 
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Expanding all the aerodynamic forces and moments of the former 
set of equations as a Taylor series in the appropriate variables, but neg¬ 
lecting the higher order terms; 


* Q. 




d/t '~y~ Sot 

M %- 'v + ••• 

and letting: 

- *r a</ 

o 

O « 9^ + & Q 

V s - Vn r AV 


cS- of + A of of * A(/ 


£ « ^ + A0 


y = V£ * ay 


B - aQ 
V = aV 


(7a) 


(7b) 


(7c) 


(7d) 


( 8 ) 


>r * # - «r 


r - 


where the delta variables denote the perturbations and a , 9 and V 

o o o 

are equal to zero, since it is assumed that the motion takes place about 
some steady-state trimmed flight condition where the initial rates are 
zero. 

Also, during the disturbed motion of the aircraft, the rotor angular 
velocity £2 is assumed constant, so that Ajx = AV/£2R . 


j t. 



Now solving the longitudinal dynamic equations (Eq. 6) in terms 
of the perturbations, the convenient substitution of Eq. 6b into Eq. 6a 
may be made: 

/rn — V/ { 

so that Eq 6a becomes: 

= - of (7~+ 1 — V!6 = o ( 

Substituting Eq. 7 and Eq. 8 into Eq. 9 and neglecting higher order terms 
the H-Force equation becomes: 

ZF X - -/*» ^ -// e -D a -we c -°f(r a -w + c<rD o )- of pV 

For the initial steady-state condition (subscript o) before any 
acceleration occurred, Eq. 4 may be written as: 

-# c -6 a -we a = c (1 

- T o -<F 0 Q, + W = c (1 

*71--x T +D 0 -^=o 


Substituting Eq. 11a and Eq. lib in Eq. 10 reduces the H-Force equation 
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0--c/JlLaU - Maju, 

0 9/^/ 0 dof 3yU 


~~ - ~ A6 — -—■ AM. - Wa6 

3& d/t / 


(12) 


2 2 

Non-dimensionalizing this equation by dividing Eq. 12 by prrR (i2R) and 


treating the drag term as: 


b 




r - . _ s _,_ 

° f> (flfi) 2 pfrfi : (Slfi) 1 2 S ^ 2" 


where: 


then: 


c '4 

*> S. 


(13a) 


(13 Jo) 


= jL (c‘ A 1 = r / i 6 ^ 2 ) = alC' 

C/b d/l ^ 0 2 / £> 9/U ? / D 


(14) 


Also introducing the aerudynamic time parameter: 

>791/ 


r ~ 


PTrfi) 2 - fuR 


and using the notation: 


if )= sL L ) 
} d(W) 


(15) 


(16) 


where time is counted in the dimensionless time ratio t/r , the H-Force 
equation in the usual aerodynamic form is: 


WW'% 
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[ C "u c pt f A c a ^ ^L r v / Cr J\ A ^ 

(17) 

+[c + C H = 0 

L » d(A<9) .! 

The thrust force equation may be developed directly from Eq. 6b using 
the same methods: 

^pn V Y - T-cZ D ■+ W - O 

(6b) 

and using Eq. 15 and Eq. 16 with Eq 8: 

* * 

dr = 

dG - d (a6) 

dof - d (Aod) ^ ^ 

the thrust force equation becomes: 


i A 


V 4 C 


i + A 


d &<>i-/U e d LG *C 


For the moment equation (Eq. 6c) non-dimensionalization is accomplished 

2 2 2 
by using pnrR (SIR) R , and defining mk Y or; 



where k is the radius of gyration and h^ is the non-dimensional 
moment of inertia parameter. 








17 


The moment equation is: 

b*. * c r: ^ - C' . 


?■ /* >? 


v * 


>' ^ jfck*-NV 


CL + C r . £r 

c*T /? 


C 


• d - C • %-<z-Ad%e = o 
’ho) di*&) * r j 


( 20 ) 


Hence, the three equations of motion in familiar aerodynamic form are 
Eqs. 17, 13 and 20. 

H-Force Equation 


' r 
'H 


A< 


- ^ c r 
c > 


^ s - dU A + [y + y <y ^ 


+ [ C v * C "e(, 6 ) -^U& = O 


d 7 ) 


Thrust Equation 

fe c °w 

Moment Equation 


it / r 






(18) 


& w c: g. 


r 




( 20 ) 


C 

* V yf 


Uo^ t 
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Explanation of the Terms of the Equations of Motion 

The three equations of motion are simultaneous, homogeneous, 
ordinary, differential equations with constant coefficients. These constant 
coefficients are expressed by stability derivatives, velocity ratios, a 
non-dimensional inertia parameter, angle measurements and drag coeffi- 

c c c 

cients. In this Report, the stability derivatives H , H , T 


^m a«d 

a a 


M <*' f 

are evaluated from data obtained by wind tunnel 


experiment. The method of evaluating the stability derivatives and drag 
terms is explained in the Appendices. 

The individual terms, expressed in non-dimensional or coefficient 
form, in the equations of motion, in effect, represent the following: 

Term 


1C H 

M 


2. C h -A a 

a 

3. C T * A ii 


4. C T • A* 

a 


represents the H-Force due to the change Ap in the volo- 
C 

city ratio where H is the rate of change of H-Force 

^ C 

coefficient with velocity ratio. For the case of H 

P 

H-Force obtained from wind tunnel data included the 
measurement of all the aerodynamic forces of the rotor 
and duct assembly acting normal to ‘he thrust axis, 
represents the H-Force due to a change A a in the angle 

Q 

of attack where H is the rate of change of H-Force 
a 

coefficient with angle of attack. 

represents the thrust force due tc a change Ap in velo- 

Q 

city ratio where T^ is the rate of change of thrust force 
coefficient with velocity ratio. 

Q 

For the case of T , thrust force obtained from the 
wind tunnel data, included the measurement of all the 
aerodynamic forces of the rotor and duct assembly acting 
along the thrust axis . 

represents the thrust force due to the change Aa in 

Q 

angle of attack where T is the rate of change of thrust 

a 
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6 - “VoS' 

2 

7/f C n 'Aff 
2 D 

8. C • 
m 

/i 


9. C • Ao- 
m 

a 


force coefficient with angle of attack, 
represents the drag force in the X direction contri¬ 
buted by the pilot due to the change ap in the velocity 
ratio. The term C ' equals --and p is the initial, 

D S R 0 

steady-state velocity ratio. 

(See Appendix I for additional explanation of terms 5,6, 
and 7.) 

represents the drag force in the Z direction, contributed by 
the pilot, due to the change Ap in the velocity ratio. The 
initial, steady-state angle of attack is defined aso^. 
represents the drag force in the Z direction contributed 
by the pilot due to the change A a in the angle of attack, 
represents the moment contribution measured about the 
quarter chord point, due to the change Ap in the velocity 

Q 

ratio where m is the stability derivative denoting the 
A* 

rate of change of pitching moment coefficient about 

the quarter chord point with velocity ratio. 

£1 

The stability derivative *m was evaluated in wind 

b 

tunnel experiments by measuring the moment created by 
the duct and rotor only. This moment was arbitrarily 
transferred to the quarter chord point for convenience 
in accomplishing additional transfers when vertical c.g. 
location was varied in order to study the effect on the 
stability derivatives . Changes in vertical c.g. loca¬ 
tion change the distances e^, x^, and a^. 
represents the moment contribution measured about the 

quarter chord point due to the change Aa in the angle of 
C 

attack where m^ is the rate of change of pitching moment 
coefficient about the quarter chord point with angle of 
attack. This derivative was evaluated as explained in 
term 8. 


i. 
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12. ^C H -A„ 

* d 


13 ■ T c h 

a 


D 

14. —p C 'Au 
R D M 


15. C W A6 


Ap represents the moment contribution about the 

C 

center of gravity due to the thrust force T • Ap 

h 

acting at the moment arm x^. 

represents the moment contribution about the center 

c 

of gravity due to the thrust force T • A a acting at the 

a 

moment arm x^. 

T 

represents the moment contribution about the center of 

Q 

gravity due to the H-Force • Ap acting at the 

moment arm z . 

H 

represents the moment contribution about the center of 

Q 

gravity due to the H-Force H -A a acting at the moment 

a 

aim z . 

H 

represents the moment contribution about the center of 

9 

gravity due to the drag force p q C d ' Ap acting parallel to 
the relative wind. The perpendicular distance from the 
center of gravity to the line of action of this drag force 
is the moment arm z^ 

represents the weight force component in the X direction 

due to a change A0 in attitude angle of the machine. For 

the component C,.,9 in the X direction, C,.,0 drops out 
W Wo 

from Eq. 12 due to the initial condition (Eq. 11a) so that: 


C 6 = C 9 + C. A0 
W W o W 


or: 


C 8 = C A6 

W W 


However, in the thrust equation, the weight term has no 
perturbation component and drops out entirely from Eq. 
6b due to Eq. lib. 


awH 'i m 'wiw min i tin rani 


T.~> nwnMM 







21 


16. C. 


d(A 0 ) 


■d(A@) represents the additional H-Force due to the linear 
velocity caused by an angular pitch rate d(A0) about the 


17. 2 H 


center of gravity. 

c 

The derivative H represents the rate of 

d(A 0 ) 

change of H-Force coefficient with pitching velocity 
d(A£). 

^d^0'" re P resent:s the moment contribution about the center 

of gravity due to the H-Force . /A ^.acting at the moment 

d (A 0 ) 

arm z . This term can be thought of as a contribution to 


18. C 


*d(A 0 ) 


19. d(A^) 


the damping in pitch of the aircraft. 
d(A0) represents the moment contribution due to the pitch 
rate d(A0). This term is the familiar pitch damping term 
and in these equations refers to the pitch damping cre¬ 
ated by the rotor and duct. Determination of this value 
is difficult but an attempt has been made to evaluate the 
rotor and duct contribution to pitch damping in Appendix II. 
represents the component along the X axis of the inertia 
term due to a linear acceleration of the aircraft. 


20. ^ d(A a) - ij. d (A 0) ■ 
o o 


•taken together, these terms represent the 


component along the Z axis of the centrifugal force 
inertia term caused by the rate of change of flight 
path angle 7 . 

21. A p + A a -taken together, these terms represent the 

component along the X axis of the centrifugal force inertia 
term caused by the rate of change of flight path angle 7 . 
The form of this term is due to the substitution performed 
in Eq. 9 where £q. 6 b: 


sKi V f = T + - W 


( 6 b) 
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wad substituted into the H-Force equation as a matter 
of convenience in solving the set of equations. If 
Eq. 6a had been solved as it stands, this inertia term 
would have been in a form similar to term 20 in the 
thrust equation. However, the centrifugal force compo¬ 
nent in the X direction, mV?#, of Eq. 6a, after the mathe¬ 
matical substitution and operations, takes the form of 
term 21: 

V L "" ^ - sin VfoC 

This indicates that this bodv axis inertia term is now 
represented by components of forces parallel to the rela¬ 
tive wind. All terms in the original equations represent 
the summation of forces and moments along and about body 
axes and, in this case, because of the aforementioned 
substitution, the body axis inertia force is represented 
by components of forces parallel to the relative wind. 
However, these still represent only the body axis inertia 
term and the equations are still consistent with the body 
axis frame of reference. Changes in variables and sub¬ 
stitutions among the set of Eq. 6 in no way change the 
validity or meaning of the original body axes equations 
of motion describing the dynamic system. 

2 

22. h y d“(A9) represents the pitching moment inertia term due to the 

distribution of the mass of the vehicle where h^ is the 

non-dimensional moment of inertia parameter measured 

2 

about the Y axis. The term d (A0) represents the per¬ 
turbation angular acceleration about the Y axis measured 
with respect to the time ratio (t/ T ) as indicated by the 
symbol (d). 
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Solution of the Equations of Motion 

In order to solve the equations of motion the form is altered 
slightly by letting: 




r 


Cl 


(21a) 


( 21 b) 




< 




C 



< 5 , ^ 


( 21 c) 


(2 Id) 


C- C m +*r C r 




^C, 


(2 i e) 



£ 

/? 



(2 If) 




<ZkO) 


= r 


^ 6 >; 


— ^ 

* <£(*e) 


(2 lg) 


where the bars over the stability derivatives indicate the overall or total 
values of the forces and moments of each specific derivative acting along 
or about the axes . 

From the above notations it is possible to see that changes in 

C C 

vertical c.g, location will affect the values of m and m because of 

M o 



changes in z TT , z , and x m . The term ,, A ^ will not change radically 
c H d T d(A6) 

since is a rather small quantity. 

Also, the moment of inertia about the Y axis L^, is held constant 
for this dynamic analysis. 

Using this notation, Eqs. 17, 18 and 20 take the form: 

H-Force Equation: 




• d. \ & @ — o 

d(a9) J 


T-Force Equation 


[%]*/* ° 


Moment Equation 


i C -J A ^ + [%^ d " *r 1 *} e ‘° <• 

The equations of motion are three simultaneous differential equa¬ 
tions with constant coefficients, the solutions of which may be obtained 


by assuming a solution in the form: 

, , & 

A/U =■ Atf, e 


aO=a6 / c 


/ 

d. (^o)- AfyJe 


2 

d (&6)= aGjA e 


where X is a real or complex constant of equal value for each variation, 
and where A^, Ao^, and A9 , are also real or complex constants. 
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Substituting Eq. 23 and Lq. 24 into £q. 22, and dividing out 
X t/j 

the common term e ' , reduces Eq 22 to three algebraic equations 

in the unknown X and the terms A^, Aq^ and A0^; 


r^ 1 


t + C^-f Q, 'O (2oa) 


[% 


(25 b) 


L /i\ 




+■ 


c 


/*n 


£*£ 4- 

i 


S 

L d(d&) 


) - hy } JAfi> = o 


( 2 ^c! 


It is a property of this type of equation that there can be non-zero 
values of the unknowns only when the determinant of the coefficients is ° 
equal to zero: 






+ £ 9 X 


C, 




<L 

y9 yu. 


Cr AA° x * 


c 


£ -J-Ayl 

d^O) Y 


= 0 


(26) 


Expanding this stability determinant; the characteristic equation 
in X is of the form: 


A\*+BX*+cf +&X+E — ° 


(27) 



Beet.- 
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where: 


/T 7 ••= Ayr/Lo 


(28a) 


B- 


% 





C = 

- 4 







</<ie) 

£ 3 ~ ^ 

if 

y > *v- 

0 = 


c 

/* 

'"’ike) 

+/*• 

. c c c 


+/ t o 

C w 

£ 


c 



-A> 

cL 

7* 


& "oS 

5 

6) 


£= 



- C r ' 

1 



( 28 b) 


(28c) 


(28d) 


(23e) 


It is convenient at this point to consider the hovering Flying 
Platform as a special case of the former development. 

The common assumption, in hover analysis, is the exclusion of 
the vertical degree of freedom since, in hovering, the horizontal motion and 

c c 

tilting do not affect the thrust. The derivatives T and T are 

[1 a 

are zero. Hence, the thrust will be essentially constant and there will 
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be no appreciable vertical motion of the machine. Principal interest 
is centered cn changes of velocity ratio and attitude angle and not on 
altitude change. This can be accomplished mathematically by altering 
the original equations of motion so that the thrust force equation is elim- 

inated entirely. The derivatives and are zero and the alpha 

te-ms in the H-Force and moment equations are dropped. 

The algebraic equations are: 

{S" / c+X)Y', -*-(C w + = o 

_ _ (29 

(?y,) A A‘ + - tyj 2 )*e, = o 

and the:determinant of coefficients is: 








= O 


Expanding this determinant yields the third order characteristic 


equation: 


CJ + D-o 

rr rr n ft 


where: 




(32a) 
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B - /? T CL ~ (TL 

c = c c - c c 

* '>*<£(*&) V* '*’*£*) 

D - C C 

# ^ -H' 


(32b) 


(32 c) 


(32d) 


The characteristic equations obtained for forward flight and 
hovering are often referred to as the stability quartic and cubic. The 
roots of these equations are the values of X that determine the final solu¬ 
tion. 

Then, for example, the response of an aircraft in pitch AG to a 
disturbance may be found by a solution of the following form: 






A6 = Ae,e *-A6 2 e /A&e + aB v c 


(33) 


Determination of time histories of the variables by actual 
solutions of this form is quite lengthy and tedious if numerous deriv¬ 
ative variations are made. 

However, analog computers can be used to solve the equations of 
motion in order to obtain the time histories associated with arbitrary 
changes in the stability derivatives. Typical time histories are shown 
in Figs. 5 and 6. 
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Investigation o f the Roots of the Char act eristic Equation 

Frequently, the important aspect in the study of a dynamic system 
is the determination of the characteristic modes of motion of the aircraft 
after a disturbance from steady flight, by an investigation of the four 
roots of the stability quartio. 

Methods of predicting the character of the motion from the roots 
or the coefficients are covered in various texts. (Ref. 5). 

It is of interest to know whether the motion is oscillatory or 
aperiodic. If oscillatory, it is necessary to determine the period and 
datnping. When aperiodic, it is important to determine the rate of 
divergence or convergence. 

This information may be obtained by examining the roots of the 
stability quartic. 

When the roots are positive real numbers, the motion is an aperi¬ 
odic divergence; negative real numbers, aperiodic convergence; and 
complex pairs, an oscillatory motion. In the case of a complex pair of 
roots, the oscillation is undamped if the real part is positive, and damped 
if the real part is negative. 

The roots of the quartic may be obtained by any of the numerous 
methods outlined in many texts. 

In general, if the roots are both real and complex as. 


= Q,±ilr 

h = e 

Ay ~ Oi 


(34) 


the period of the oscillatory pair may be found by; 



Zjrjr 

lr 




[seconds] 


(35) 
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and the damping may be found as the time (T) required for the transient 
response to grow to double amplitude or to decay or damp to half 
amplitude: 

——ft- 7" =• —■ rS£CC/?&] 

a CL L J (36) 

Taking the center of gravity at the original location, (.42c above 
the top of the duct) the roots of the quartic for the different ducts at 
different speeds are computed. The damping and period are shown for 
each of the roots. See Tables I and II. 

I nvestigation of the Coefficients of the Characteristic Equation 

Considerable information concerning the character of the motion 
of the aircraft can be obtained from the coefficients of the stability quartic 
without solving for the roots. 

If all the coefficients of: 

+• = O (37) 

are positive, there can be no positive real root and there is no possibility 
of a pure divergence. If one of the coefficients changes sign (becomes 
negative) there cai) be either an increasing oscillation or a pure diver¬ 
gence in one of the modes . 

The C coefficient is of special importance since it represents 
approximately the maneuver margin of the aircraft (Ref. 6). If the C 
coefficient changes sign (become^ negative) it is an indication of a 
serious instability. There will probably be a rapid divergence or unstable 
oscillation. 
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The E coefficient, where: 

r- = SI l~ iCr _ d£m dCr 

* L <9A dot dot d/U,^ (38) 

is of importance for static stability. The static stability is positive if E 
is greater than zero. If E is equal to zero, one of the roots is zero and 
one of the modes can continue unchanged indefinitely. Usually E is the 
first coefficient to change sign, and therefore the first indicator of a 
divergence or positive real root. In some instances, it is possible for C 
k to change sign before E, in which case there will be a dynamic instability 

r 

occurring with a positive static stability. 

Regardless of the signs of the coefficients, it is possible to deter¬ 
mine the number of roots having positive real parts by using Routh's 
criterion. This method requires only simple calculations utilizing the 
coefficients of the stability quartic in setting up the familiar Routh's 
array. The number of roots having positive real parts is equal to the 
number of sign changes occurring in the first column of the array. 

Also, from Routh's array, it is possible to obtain the combination 
of coefficients known as Routh's discriminant: 

R = BCD -/ID 1 -B* E (39) 

where A, B, G, D and E are the coefficients of the quartic. If all the 
coefficients are positive , it is possible to make the followino statements 
concerning the stability of the system. 

If Routh's discriminant is positive, there is no possibility of the 
real part of any complex pair being positive. There will be no undamped 
oscillations. 


• -T-unrnrriwic 
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If Rouih's discriminant is negative, there is a complex pair of 
roots with a positive real part. There is an undamped oscillation. 

If Routh's discriminant is equal to zero, there is a ne>jiraity'damped 
oscillation implying a complex pair of roots with a zero real part. 

The criterion, (R = 0), then indicate? a boundary between stable 
and unstable oscillations. 

This condition, (R = 0), developed from the coefficients of the sta¬ 
bility quartic, may be investigated further by,assuming a solution of the 
form for neutral oscillations (A = iu). Substituting this solution with the 
zero real part into the characteristic equation gives: 

* CJ%DJ. +£ = o 

At cv Y'Si 3 ^ t- Ci 2 w 2 + Dtcu /- S - c ( 4 °) 

&u> 3 i - Ccu* +&tui i- F = o 


Forming two equations of the real and imaginary parts: 

/?u> * — C u> * ■r- - o 

- £oj* + Du - o 











Substituting Eq. 43 in Eq. 41 results in: 


de. _ . , „ 


£C£> ~/\D 2 -B 2 £ = 0-/9 


which is the form of Routh's discriminant. 

With this analysis, one can obtain an appreciation of the require¬ 
ments imposed on the coefficients A, B, C, D, E, so that the solution 
R = 0 is a stability boundary between stable and unstable oscillations. 

Since a solution of the form for neutral oscillations is assumed 
(X = iu), inspection of Eq.-43 reveals that Ji will be a real number, e.g. , 
u = a, if D and B are both positive, and the solution X = ia will 
indicate a neutral oscillation. However, if D or B changes sign so 
that w = ia in Eq. 43, then the solution will be X = -a, or a real root. 
Therefore, when the D coefficient changes sign, (becomes negative), 
the solution of Routh's discriminant equal to zero no longer defines an 
oscillation boundary. 


Boundaries of the Stability Quartic 

In_the analog computer program, it was found that the derivatives 

c c 

m and m had a principal effect on the mode of motion of the Flying 

hi a 

Platform. 

The.Figs. 10A, 11A, 13Aand 14A entitled "Boundaries of the 

C C 

Stability Quartic" are drawn with m and m as the ordinate and 

ja a 

abscissa. 



i *' V •" 
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Included bn this diagram are various boundaries of the stability 

. Q Q 

quartic plotted as functions of m and m . Values of the total stability 

yX Q 

derivatives were used as calculated from the wind tunnel data for each 
speed and duct configuration. The theoretical value of pitch damping 
was used and kept constant for these diagrams. 

The coefficients oi the stability quartic, G, D and E are plotted 

c c 

cis functions of m and m and appear as the straight lines labeled 

fi a 

C = 0, D = 0 and E = 0. The lines so labeled represent the boundaries ajt 
which each coefficient changes sign from positive to negative for the 

c c 

appropriate values of m and m . 

IL a 

The coefficients are negative on the hatched side of the line. For 
example, near the origin in the right half plane the vertical line labeled 
C = 0 is negative on the right aide of the line. 

The meaning of the traversing of these lines (coefficients changing 
sign from positive to negative) has been discussed previously and, for 
this section, can be summarized as: 

C< 0 - Serious instability associated with negative 
maneuver margin. 

D < 0 - The solution R - 0 no longer has meaning as 
an oscillation boundary. 

E"< 0 -Instability associated with a negative static 

margin. Aircraft statically unstable. Dynamic 
instability present in the form of a divergence 
due to static instability. 

c c 

The coefficients A and B, not functions of m and m and 
generally positive, do not appear on these diagrams. 

The line R = 0, representing Routh's discriminant, is plotted 
to denote the boundary between stable and unstable oscillations. 
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C 0 

This line represents the loci of the values of m and m for 

a 

neutral oscillator/ motion. 

On die left Side of the line R is positive and the oscillations will 
be stable . On the right or hatched side R is negative and the oscillations 
will be unstable. The line R = 0 does not exist past the line 3=0, since 
its meaning as an oscillation boundary is ended whenever the D coefficient 
becomes negative. 

In one case, (Fig. 13 A), the line R = 0 has essentially three 
branches, all representing a neutral oscillation boundary, and all located 
in the region where the D coefficient is positive. In the top half plane, 
the two branches of R = 0 form a parabola ineide of which there is an 
unstable oscillation (at least two of the roots must foric a complex 
pair with a positive real part). In the bottom half plane, only one 
branch of the parabola appears since the other is nonexistent when it 
crosses the D = 0 line into the region where the D coefficient is negative. 
Note that the line E = 0 passes through the middle of the R = 0 parabola in 
the top half plane. 

c c 

Whenever the values of m and m locate a point on the negative 

H a 

side of the E = 0 line, at least one of the roots of the quartic must be 
positive end real, and an instability exists which takes the form of a 
divergence. 

Then, in the region between the E = 0 line and the right hand branch 
of the R = 0 line in the first quadrant (Fig. 13A), at least three of the roots 
of ttfe quartic are accounted for; a complex pair with a positive real part 
and a positive real root. The other root in this region is a convergence 
as seen in the root It cus plot Fig. 13-1. 

The other two branches of the R = 0 line located in the second and 
third quadrants represent oscillation boundaries , so that in the region 
between them, R is greater than zero and the entire area represents a 
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completely stable region of only convergences and stable oscillations. 

In Fig. 14A, the R = 0 line is terminated when it reaches the E = 0 

Q 

line since it passes into a region where m^ and E are always negative, 
so that the predominating mode of motion is a divergence that strongly 
masks all other modes of motion. 

The line labeled RLBP represents a locus of points obtained from 
root locus diagrams. The letters RLBP are an abbreviation for "Root Locus 
Break Points". The break point in this case represents the point in the 
right half plane of the root locus diagram where a complex pair of roots 
joins the positive real axis, thereby terminating the unstable oscillation 
and becoming two aperiodic divergences. 

In other words, moving to the right from the stable region, across the 
R = 0 line , and passing into a region of unstable oscillations, at least 
two of the roots of the quartic are a complex pair with positive real 
parts (Figs. 10A, 11 A, and 14A). Continuing to move to the right.into 
the first quadrant and crossing the RLBP line, this Complex pair becomes 
two real positive roots, ending the unstable oscillation and now repre¬ 
senting two divergent roots. To the immediate right or hatched side of 
the RLBP line, at 14a st two of the roots of the quartic are positive real 
roots. 

The line RLBP does not appear in Fig. 13A, because the parabola 
R = 0 already defines a large region where unstable oscillations start 
and stop. 

Q Q 

However, for extremely large values of m and m beyond the 

/i a 

values of this diagram (Fig. 13A), there is a region where the complex 
roots become positive real roots for a time. (See root locus Fig. 13-I-e, 

c c 

points 180 and 181.) The value of m is equal to 0.006 and m is 

O' /i 

approximately equal to 0.1 at these points, which ire beyond the pertinent 
C C 

rang# of m and m for this particular case. 
a 


vet 











37 


The various center of gravity locations have been plotted on these 

c c 

diagrams according to their respective values of m and m as 

/L a 

listed in Table III . The broken line between these points represents the 
approximate line,- of travel for the c.g. at intermediate values. 

Explanation of Stability Diagrams 

The Figs. 10B, 11B, 13B and i4B entitled “Stability Diagrams" 
are a greatly Simplified form of the previous boundary diagrams. The 
only portions of the boundaries that have been retained are those which 
divide the diagram into the pertinent regions of stability. 

The unshaded area represents a region of complete dynamic sta¬ 
bility . Only roots representing convergences and stable oscillations 
exist in this region. 

The shaded areas represent regions of dynamic instability. Roots 
existing in this region represent divergences and/or unstable oscillations. 

There are always four roots to contend with in every region. In 
the stable region all four roots must be stable roots. In the unstable 
region, a minimum of one unstable root must always exist while the others 
may be stable. 

Actually, the unstable roots are the predominant modes of motion In 
the unstable regions and mask the stable modes present. 

One shaded unstable region is labeled "unstable oscillations" 
and indicates that at least two of the four roots in this region must be 
a complex pair with a positive real part. The other shaded area indicates 
the unstable region where a minimum o f one unstable root must exist. 

Some of the interesting modes of motions which prevail in the 
various regions are labeled in these diagrams. 






Again, the various center of gravity locations have been plotted 

c c 

according to their respective values of m and m as listed in Table 

a 

III. The broken line represents the approximate line of travel of the center 
of gravity at intermediate values. 

C 

For the case of hovering, u. and m equal zero. The diagram 

C ° a 

is plotted only for m and appears in the form shown in Fig. 8. 

The study of these diagrams, together with the associated root locus 
diagrams, will present quite clearly and concisely the precise modes of 
motion to be expected at different speeds for different center of gravity 
locations with the different duct configurations. 

f 

Explanation of the Root Locus Diagrams 

The stability diagrams were prepared in order to show concisely 

c c 

and clearly the role of m , m , vertical center of gravity travel, speed 

/i. a 

changes, and different duct designs, in affecting the modes of motion 
of the Flying Platform. 

The diagrams exhibit rathter precisely the mode of motion associated 
with locations in the various regions. 

Although the determination of the stability boundaries is quite 
lengthy, they may be calculated by various methods. 

In the formulation of the stability diagrams for this Report, use 
was made of the root locus method of displaying dynamic responses. 

With this method, it is possible to exhibit the properties of a linear 
system by plotting the loci of the roots of the characteristic equation in 
which the values of selected derivatives have been varied. 

A complete treatment on the utilization of this method is avail¬ 
able in texts by Evans (Ref. 7) or Savant (Ref. 8). 
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Since primary interest is centered on varying the derivatives m 

C 1 

and m , the characteristic equation may be written in the form: 
a 



where the quantities enclosed in the parenthesis are equal to the coeffi¬ 
cients C, D and E. From this expression it is possible to calculate the 
polas»nd the zeros. Since vertical center of gravity travel seems to have 

c 

the primary effect of changing m , the root locus diagrams were plotted 
C ^ 

so that m varied from zero to positive or negative infinity. The poles 

M C C 

were calculated for specific values of m with m equal to zero. 

a il 

Solving the characteristic equation in this manner gives the four roots of 
fhequartic p , p , p , p . The zeros of the equation are the increments 

Q J. Z 0 T 

of m , namely: 


C 3 

d 3 ( c » 


4 


= c 


tv 



( 48 ) 


Q 

Solving the above, the zeros for m are obtained and denoted by z. 

C ^ 

Similarly, the zeros for m denoted by Z, may be obtained from: 

a 


c , - 

D .> = - C r C„ -/*.& 


( 49 ) 


£ 
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Using the former expressions, the root locus equation may be writ- 

c c 

ten for a specific value oi m with m varying from zero to positive 

a p. 

or negative infinity as: 

/'o^y feoYCw + E f ) = 0 ( 50 ) 

c 

where the p's represent the poles and the z's the zeros of m . 

C P Q 

For the case where the poles (?) are calculated for m = m = 0, 
C M a 

and the zeros (Z) for are utilized, the expression takes the form: 

4 - c m Z)U-K) - 0 < 51 > 


c c 

in which m^ varies from positive to negative Infinity while m^ equals 
zero. 

A group of root locus diagrams are plotted so that essentially all 
C C 

values of m and m are covered for a particular duct at a selected 

pa 

speed. 

Associated with each stability diagram there are five root locus 
sets. The first set investigates the modes of motion encountered when 
moving along the ordinate or abscissa of the stability diagram. The 
other four sets, labeled I, II, III and IV, investigate the modes of motion 
present in the four quadrants of the stability diagram. The value of 

Q 

m^ (Eq. 50) in these diagrams varies from zero to positive or negative 
infinity. 

The root locus sets are plotted for both ducts at pi equal to 0.05 
and 0.10. 
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Included In each figure used to investigate the modes of motion in 
each quadrant (I,_II, III, IV) , one root locu3 is plotted with the zeros of 

c c 

m instead of m . This plot is divided from the rest by a double line 
a 11 

border as in Fig. 10-I-f. In each of these particular double border 

c c 

diagrams, m is equal to zero and in varies from zero to infinity. This 
hi a 

then shows the movement of the initial positions of the poles for each quad- 

C 

rant as the value of m is increased from zero to infinity. In this way it 

* c 

is possible to sketch the root loci for all higher values of m beyond 

a C 

those calculated for each quadrant. In this respect, the values of m 

a 

in each quadrant are usually large enough so that the basic appearance 
of the root loci will not change for values beyond the highest shown in 
each set. 

All root loci drawn in this report were accurately calculated and 
plotted to the scale ahown in each figure. 

The root locus set for ^ = 0 is also plotted in Fig. 7. For both 

Q 

ducts in the hovering case m equals zero* The stability cubics are 

a 

written as: 




(31) 


or: 






O 


(52) 


The root locus equation, after dropping a negligible "zero" due to 

Q 

the small value of H, , * , is: 

d (AQ) 




( 53 ) 
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where C^. Is known and mis varied from zero to positive or negative 

Q 

infinity. The poies are calculated for m equal to zero. The diagrams 
for the hovering case have three poles and no zeros. 

The root locus method may also be utilized conveniently for the 
investigation of the effects of artificial stabilization on the transient re-- 
spCnseof the Flying Platform. Artificial stabilization may be provided 
by a rate sensitive device such as the gyro-bar with aerodynamic or 
mechanical damping. 

The addition of this type of device to the dynamic system may be 

accomplished mathematically by ir eluding a rate sensitive transfer function 
;> Kd 

with a phase lag of the type y—-g- — to the pitch damping term in the 
moment equation. 

For hovering, the determinant of coefficients in algebraic form 
becomes: 



Expanding this determinant: 


0 = A y /l + (by ''•V +JZL (x+C# ) (55) 


7* &6>) 

and arranging the poles and zeros as: 


o = ['^ t -0D> v /?W'vC)^-C C, a + C„c 

L JL Y a Y ^ 


/ (56) 




43 


where T is the time constant for the phase lag , the root locus equation for 
the hovering case may be written as: 

4 - ) = o ( 57 ) 


where the gain of the feedback loop, K, of the stabilizing device varies 
from zero to infinity. 

The determinant of coefficients for the addition of artificial stabili¬ 
zation in forward flight is: 



c 

*eC 

c +c . 
* 

•A 


l -p& 

<5, Y*° J ' 



= 0 

( 58 ) 

c 

d 


ArA 

// TA 



and expanding the determinant: 






/a 




( 59 ) 


where the first five terms represent the original characteristic equation 
and T is again the time constant for the phase lag. Separating the poles 
and zeros in root locus form gives: 


^ ^ },)U- &) = O 


> 

i 
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where K, the gain, varies from zero to infinity. 

The root locus diagrams for the artificially stabilized machine in 
hovering and high forward speed are plot f ed in Figs, 9, 12 and 15. 

These diagrams , together with all other pertinent root locus 
diagrams, are fully explained in the DISCUSSION. 


Effect of Vertical Center of Gravity Movement oh the Overall Stabil^ v 
Derivatives 

The overall stability derivatives are expressed by £q. 21 as: 

(2 If) 

(2 le) 

The distances x^, z^ and z Q represent the moment arms from 

the total center of gravity to the respective forces, ^m and ^m 

/l a q 

lepresent the stability derivatives about the quarter chord, whereas, m 

C ^ 

and m are functions of the distances x_, z„ and z^, When the total 
q- T H D 

center of gravity of the aircraft is moved vertically these distances will 
vary so that considerable change* in the overall value of the derivatives, 

c 

especially m , can be obtained. 

M 

In the original c.g. location, the total center of gravity is 0.42 
of the chord length (0.42c) above the top edge of the duct. For the case 
of the long chord duct, the total center of gravity of the aircraft is located 
approximately nine inches above the top edge of the duct. Similarly, the 
center of gravity of the short chord duct is also assumed to be at the 0.42c 
location so that it is approximately five inches above the top of the duct. 
These positions are designated as the original center of gravity positions 
from which the vertical variations of c.g. are measured. The pilot's 






A 

* ^ 




4 C 


c 
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A 


£ 








sta id-on platform always remains at this original position, i.e., 0.42c 
above the top of the long and short chord duct. 

The actual scale profiles of the long and short chord duct models 
are shown-in Fig. 2. 

c c 

The values of the stability derivatives m and m for the different 

li a 

center of gravity locations shown in Fig. ]6 are presented in Table III. 





Pilot's 

"Platform 



c.g. 0.5 ft. original 

down froffl location 

original of c.g. 


c.g. 0.5 ft. up 
from original 
location 


c.g. 1.0 ft. up 
from original 
location 


Fig. 16 CENTER OF C 7ITY LOCATIONS 
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DISCUSSION 

General 

1. Similarity of Stability Derivatives 

From the information contained in Table III, concerning the over¬ 
all value of the stability derivatives, it is apparent that there is some sim¬ 
ilarity or overlapping of values for the two ducts for different speeds at 

Q 

different c.g, locations. For example, the values of m^ for the long 
chord duct at p equal to 0.05 with the c.g. raised 0.5 feet is equal to 

Q 

0.0136 while m^ for the short chord duct at n equal to 0.05 with the c.g. 

at the original location is approximately the same or equal to 0.0121. 

A small adjustment in c.g. location of one duct configuration would make the 

derivative equal that of the other configuration. Similarly, the value of 
C 

m^ equal co 0.0032 also matches up for both ducts when /i equals 0.10 
at the highest c.g. location. 

Therefore, at a specific speed in a certain speed range, similar 
values of either stability derivative may be obtained for both duct con¬ 
figurations by adjusting the vertical center of gravity location. 

2. Attitude Angle for Steady Forward Flight 

As seen in Fig. 18, the attitude angle for steady-state equilibrium 
in forward flight is plotted versus /*. For convenience, Fig. 17 is plotted 
directly above Fig. 18 so that conversions may be easily accomplished 
between and V in feet per second or miles per hour. 

From Fig. 18, it is evident that the long chord duct requires larger 
tilt angles than the short chord duct. This can be explained by the fact that 
the parasite and momentum drag of the long chord vehicle is probably 
greater than that of the short chord duct. 

For a long chord duct, the flow tends to align itself more axially 
than for a shorter length duct, and therefore, the momentum drag is greater. 

In order to overcome the larger total drag, the thrust vector of the 
vehicle must be tilted sufficiently so that the component of thrust in the 
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horizontal direction exactly cancels the total drag in steady-state flight. 
Although the short chord duct yields the more desirable smaller tilt angles 
other characteristics such as the performance of the long chord duct 
would hiso be considered in any design (Ref. 2). 


Static Stability 

1. The E Coefficient 

As mentioned previously in the DEVELOPMENT, the fundamental 
mathematical condition for static stability is that the E coefficient (the 
term independent of X in the stability quartic) be positive. 

The E coefficient: 


^ ~ C* S') 






(38) 


is a function of the two static stability partial derivatives m and 

C ii 

m . 

° Q 

The term v m is referred to as the static velocity stability and 

C * 

m the static angle of attack stability. 

O' 

An aircraft is statically stable with velocity,if, when acted upon by 
a pure velocity disturbance, it creates a moment about the center of gravity 
tending to restore the aircraft to the initial velocity. Similarly, an air¬ 
craft is statically stable with angle of attack, if, when acted upon by a 
pure angle of attack disturbance, it creates a moment about the center of 
gravity tending to restore the aircraft to the initial angle of attack. The 

Q 

criterion for static velocity stability requires that m be greater than zero 

C ^ 

while static angle of attack stability requires that m be less than zero. 

a 


i 
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Generally the Flying Platform has too much velocity stability in 
hovering and forward flight. However the velocity stability does 
decrease as flight speed increases and, depending on c.g. location, may 
even become unstable by changing sign. 

Q 

Wi h respect to the angle of attack stability, in hovering m 

a 

is equal to zero and the Flyjng Platform is neutrally stable with angle of 

c c 

attack. In forward flight, m has the wrong sign, that is, m increases 

a a 

positively as flight speed increases so that the platform is statically unstable 

with angle of attack. Neither of these criteria, taken by itself, is 

adequate to predict the static stability of the Flying Platform. The commonly 

C 

used static stability criterion of m < 0 for airplanes, is not sufficient for 

C a 

this case since it assumes m is equal to zero. 

For an aircraft operating in a flight regime or configuration where 

Q 

m is not equal to zero, the general criterion for static stability must 
^ G G 

include the effects of both m and m as exemplified by the E coefficient. 

fL a 

Static instability results in a dynamic divergence whenever this coefficient 


becomes negative. . In the speed regime considered for this Flying Platform, 

c c 

all the terms in the E coefficient are always positive except m^ and T . 

For the long chord duct, the term ^m^ is generally positive and large, 
but can be reduced and even made negative in certain forward flight speeds 

Q 

by raising the center of gravity. The term T becomes negative rather quickly 

C M Q Q 

and since m^ is always positive the quantity - T^ m^ is usually a 

positive increment to the E coefficient. From this it is evident that the E 


coefficient of the long chord duct, in the region of speeds considered, will 


probably_not change sign unless the center of gravity is raised sufficiently 

Q 

so that m can become a negative value. 
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Also, It is interesting to note for this case that the larger the posl- 

C c c 

tive (unstable) value of m the laraer the value of tne quantity - T . 

a 1 H a 

This becomes a larger positive increment to the E coefficient and therefore 

c 

enhances the static stability. However, large positive m has a simul- 

a 

taneous adverse effect on the other coefficients of the stability quartic. 

C 

For example, for large m , positive static stability may exist while in 

a 

all probability the C coefficient will have already changed sign, indicating 
a serious instability. See Table I. 

c 

For the short chord duct, the values of m are considerably 
smaller and become negative more quickly than the long chord duct as 

Q 

forward flight speed is increased. Also, remains positive for the 
speeds considered. The first term in the E coefficient becomes a negative 

Q 

increment as soon as m becomes negative. The second term is always 

^ C 

positive so that E changes sign as sopn as the value of m has reduced 

c c c c ^ c 

to a small positive number where T m > m T . Whenever m 

IL a il u ju 

becomes negative, E has already become negative for this case. 

Although it is desirable to have a positive static stability, instab¬ 
ilities due to E changing sign are associated with velocity changes and 
are probably not as serious as the instabilities associated with angle of 
attack changes which are indicated by other coefficients of the stability 
quartid changing sign. 

2. The Moment to Trim Versus Velocity Diagram 

The static-trim curve, in which the moment produced by the 
machine in trim equilibrium is plotted versus forward speed, can also 
be used to investigate the static stability of the Flying Platform. 

In Fig. 19, the M-V curves are plotted for both ducts for the 
various center of gravity locations by using the moment equation: 


Lrv _= H c/ + T- x r ■ H - ° 


C.* 


(4c) 
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where T • x^, is the control input of the pilot, and the moment produced 
by the vehicle including the drag of the pilot is: 


/V = Z 1 /^ - + & ■ z-o 


Inspection of the curves of this moment plotted versus velocity 

indicates that the value of the moment required to trimxan be reduced 

by raising the center of gravity position. 

The static stability, or slope of the M-V trim curve, is the total 

derivative • I n all cases presented, the slope starts as positive, 

oC v 

reaches a zero value, and then becomes negative. The criterion for 

dM 

static stability is > 0. 

r 

Also plotted on the M-V diagram is the horizontal line drawn 
at 200 foot-pounds representing the approximate maximum control moment 
attainable by the pilot for trim. Presuming that a typical 200-pound pilot, 
with his c.g. located at 3^ feet, can lean forward 20 degrees relative to 
the platform, the maximum control moment he can apply would be approxi¬ 


mately: 


/i = (200)(3. 33)(s<» 20°) 


M P => 200 ft. ~ lb-s. 

By convention, the distance x^ is taken as minus when the c.g. 
is moved forward in forward flight. The scale for plotted opposite M„ 
then, in effect, represents control position for trim. 

In forward flight, the maximum forward displacement is equal 
to approximately -0.36 feet for 200 foot-pounds. The intersections of 
the curves and the horizontal line, (denoted by circles), represent the 
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maximum forward speed for which the pilot is able to trim the moment 
produced by the Platform. 

The short vertical dashes on each curve represent the points where 

9quals zero and beyond which becomes negative. These 

points are exactly the sppeds at which the E coefficient changes sign and 
becomes negative. Actually the speeds beyond the vertical dash location 
^ z O^ f represent flight in a speed regime where the Platform is 
statically unstable. 

For a statically stable aircraft, the control position for trim always 
moves forward (piiot must be positioned further forward to trim the Platform) 
as higher trim flight speeds are achieved. 

Notice on the diagram that the control position for trim , starts 
to reduce beyond these vertical dash points. In this region, the control 
position Xj, for trim (pilot position) is located at successively smaller 
negative values for greater trim speeds after having reached the maximum 
negative value at M ai 'm-,<o 


equal to zero. This should not be thought of 


as a region of so called "control reversal", in the normal sense of the 
words, even though the control moves aft to trim each successively higher 
flight speed. Usually the term "control reversal" implies quite a different 
meaning with different effects and is most frequently associated with the 
reversed control effectiveness due to aeroelastic ohencmena., when the 
control derivative, , changes sign ( 6 is the control deflection). 

O Oq ^ 

With this meaning, for example, whenever a pilot attempts to pull the nose 
of an airplane up in order to reduce speed, (control moved aft), the aircraft 
actually noses down and tends to gain speed. Fortunately, this is not what 
occurs for the Flying Platform if it is operated in a statically unstable 
configuration, or, for that matter, *o any airplane which happens to be 
operating in a statically unstable configuration. In order to go faster or 
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move the nose down, it is still initially necessary to move the control 
forward (lean forward) to obtain the desired increase in speed. What is 
different is that once the desired higher speed is achieved, and the vehi¬ 
cle is again in trim, the pilot will find the control positioned forward less 
than it was for the slower trim position (pilot leaning forward less than 
previous trim position). In order to accelerate an aircraft in a giver, 
direction, the control will always be initially moved in the same sense 
regardless of whether the aircraft is statically stable or unstable. 

To slow down, the pilot initially leans back, but when the aircraft 
is again trimmed at the desired slower flight speed, in the region where 

~cL V iS ne9atlve ' he w111 flnd he ls leaning further forward than before. 
As seen by the curves in the M-V diagram, raising the c.g. does 
reduce the control moment required for trim but, unfortunately, also moves 
this point ( E = 0) to successively lower forward speeds, thereby 

reducing the statically stable flight envelope. 

The asterisks on the diagram indicate the maximum points where 
d / 13 equai t0 zer0 - for each duct, commensurate with the 200 

foot-pound control moment. The difference between the two radically 

different duct designs appears to be on the order of only six feet per 
second: 

Short Chord Duct » 38 feet per second 

Long Chord Duct ~ 32 feet per second 

Other flying qualities of a machine similar in nature to the 
subject vehicle are discussed in Ref. 3. 

Dynamic Stability 

The fundamental requirement for stable dynamic response is that 

all the coefficients of the stability quartic and Routh's discriminant be 
positive. 


, Rfi 
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If ail the coefficients are initially positive, the first indication of 
an instaoility will normally be a divergence when E changes sign, or an 
unstable oscillation when Routh's discriminant becomes negative. 

The roots of the characteristic equations evaluated for the original 
c.g. location of each duct at equal to 0.05 and 0.10 are shown in Tables 
I and II. 

It should be remembered that the following analysis investigates 
the control-fixed or open-loop response of the Flying Platform configurations. 

The open-loop response is the resulting mode of motion of the 
Flying Platform when the pilot remains rigid and applies no control inputs 
throughout the transient response. The closed-loop response is the 
resulting mode of motion of the Platform when the pilot acts as a feedback 
by providing controlling inputs to the transient response. Although the 
Flying Platform is unstable in the open-loop it will be stable closed- 
loop if the pilot can correctly control the machine as he desires. 

The open-loop response of the Flying Platform is investigated to 
determine the effects on the dynamic stability of vertical center of gravity 
location, duct design changes and artificial stabilization. 

In this Report this may be thoroughly accomplished by referring 
to the various figures of the stability diagrams, boundaries of the quartic 
and the associated root locus plots. 

On the figures labeled "Boundaries of the Stability Quartic" and 

c c 

"Stability Diagram", the symbols represent the values of m and m 

M a 

for each particular center of gravity location. The dashed line connecting 
these circles represents the movement of these points on the diagram, 
as the center of gravity location is varied. 
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For the different cases: 

1. Long and Short Chord Duct Hovering ^ = V = 0 

The dynamic analysis of the hovering case of the subject Flying 

Platform can best be obtained by referring to the root locus diagram Fig. 7 
a and b. 

The case is somewhat simplified by the fact that we deal with a 
cubic equation and is equal to zero. 

For the hovering long chord duct (Fig. 7b), when Is equal to 
zero, the three poles represent two convergences (point 80 and point 8» 
and a aero root (point 82)As the gain of C m Is Increased positively, 

the pole at point 80 moves toward point 86 and continues as a convergence 
for all values of positive m 
C M 

Also, as m^ Is Increased positively, the two poles at point 81 
and point 82 move toward each other on the negative real axis until they 
become equal roots at the break point feint 83). Further increases in 
cause the equal roots to break away from the negative real axis and become 
a complex pair, indicating a stable oscillation of Increasing frequency 
and decreasing stability. At point 84, C m has increased to the value at 
which the oscillation becomes neutrally damped. Increasing further 
makes the oscillation unstable (point 85)_. The oscillation continues to 
Increase In Instability and frequency as ^ is varied to positive Infinity. 

The short chord duct In hovering (Fig. 7a) has a similar analysis 
and is numbered accordingly, the only difference being the respective 
gains for the particular modes ot motion. 

In order to complete the set of stability diagrams, Fig. 8 is 
presented showing the various gains of d m for the different center of 
gravity locations together with the various modes of motion. This 

diagram Is considerably different from the other stability diagrams since 
only the values of m need to be shown. 
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For the long chord duct m^ changes approximately 0.0028 per 

inch of c.g. movement. Stable oscillations start when ^m is equal to 

C ^ 

0.0003, which is the value of m at the root locus break point (point 83). 

M c 

Unstable oscillations start when m is equal to 0.0056, which is the 
C ^ 

value of m at point 84 on the root locus diagram. All the center of 

M 

gravity locations are well inside the unstable region. 

Q 

For the short chord duct, m changes approximately 0.001 per 

M Q 

inch of c.g. movement Stable oscillations start when m is equal to 

C M 

0.000C75, which is the value of m at the root locus break point (point 93) 

C ^ 

Unstable oscillations start when m is equal to 0.0024, the value of 

C M 

m^at point 94 on the root locus diagram. Again, all the c g. locations 

are well inside the unstable region. 

For the specific case of the long chord duct in hovering, the roots 
of the characteristic equation (Table I) for the c.g. in the original locatio n 
indicate a convergent mode and an unstable oscillation. The period of 
this oscillation is approximately four seconds with a time to double 
amplitude of about one second. 

For the short chord duct in hovering, the roots of the characteristic 
equation (Table II) for the c.g. in the original location indicate a conve*r- 
gent mode and an unstable oscillation. The period of this oscillation is 
approximately six seconds with a time to double amplitude of 2.23 seconds. 

Both ducts are rather unstable with regard to the open-loop response, 
however, the short chord duct is noticeably the better of the two. 

As mentioned in the DEVELOPMENT, it is interesting to investigate 
the effects of artificial stabilization on the transient response of these 
vehicles. 

The effect of such devices on the behavioi of the system is shown 
by root locus plots for several different cases. 
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eS 

2 . Long and Short Chord Duct Artificial Stabilization Added 

Hovering [X = 0 

For the short chord duct ( c.g. in the original location ) when K 
varies from zero to positive infinity reference is made to Fig. 9a. 

The artificial stabilization adds two zeros to the hovering root locus 
diagram, one at the origin and the other close by in the left, half plane. 
Arbitrarily selecting T equal to 3.333, an additional pole is located at -0.3. 

When K is equal to zero , the poles located at point 100 indicate an 
unstable, lightly damped oscillation, while the poles at point 101 and point 
102 represent two convergences. As K is increased, the poles (points 101 and 
102) come together on the real axis at the break point 103 and then bre^k 
away from the real axis to become a stable oscillation of increasing 
frequency. Meanwhile the unstable oscillation becomes stable with no 
appreciable change in frequency (point 104). However, when K is increased 
further, the oscillation remains stable but decreases in frequency (point 105) 
until the complex pair become real roots which have limiting values at the 
zeros on the real axis (point 106). 

If we make T larger, for example equal to 50, so that the pole is 
located between the zeros on the real axis, the root locus appears as in 
Fig. 9b. 

When K is equal to zero, the poles (point 112 and point 109) rep¬ 
resent convergences and the complex pair (point 110) an unstable oscillation. 
As K is increased, the poles representing the convergences (point 112 and 
point 109) move toward and terminate at the zeros (point 108 and point 107). 
The pole at point 109 becomes a zero root when K reaches infinity, whereas 
the other pole (112) always represents a convergence. 

The unstable oscillation (point 110) increases in frequency and 
becomes neutrally stable at point 111. As K is increased, the oscillation 
remains stable but the frequency continues to increase. 

As suspected, depending on the value of the time lag, the addition 
of a lagged rate feedback device is able to stabilize the hovering platform 
to various degrees. 



JC ci 

For pure rate feedback, when T in is selected equal to zero (T = 

g = - = 0), the extra pole is located at negative infinity (a = -oo). 

As K is increased positively, the poles representing the unstable oscillation 

(point 100) move directly to the break point i06 in the left half plane instead 

of looping around (points 104 and 105) as in Fig. 9a. For this case the pole 

representing the convergence at point 102 stays in the left half plane and 

moves to the left as JC is increased. 

For pure attitude feedback, when T is selected as infinitely large 

(t = ; = s * oc), the extra pole (point 109) is located at the origin (a = 0) 
a u 

and cancels the zero (point 107) located there. The root locus for this case 
appears very similar to Fig. 9b. 

Thus, the hovering Flying Platform may also be stabilized by either pure 
rate or attitude feedback or a combination of both. 

The analysis of adding artificial stabilization to the hovering long 
chord duct is qualitatively the same as that of the short chord duct. 

3. Long Chord Duct /x = 0.05 V = 16 mph. 

Referring to the boundaries of the quartic (Fig. 10A); it is seen 

Q 

that center of gravity movement has a primary effect on the value of m 

C ** 

but only changes m slightly. Actually, on this diagram, each inch 

O' - 

C 

of c.g. movement accounts for approximately 0.0031 change in m^ and 

0.000014 change in . It is evident that all c.g. locations are 

a 

on the right side of the C = 0 line and therefore the C coefficient is always 
negative. If the c.g. is raised more than six inches above the original 
position, the D coefficient changes sign (point 1), and if raised by 12 
inches (point 2), the E coefficient changes sign. 

Referring to the stability diagram (Fig. 10B); it is noted that most 
of the c.g. locations fall within the region of unstable oscillations. If 
the c.g. is raised about nine inches (point 3) the unstable oscillations 
cease and divergences occur. Continuing to raise the c.g. (point 4) causes 
E to change sign so that even the static stability is now negative. 
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A better appreciation of the modes of motion can be obtained from 

c 

the associated root locus plots. Note that the value of m on the 

a 

stability diagram is approximately equal to 0.002 and corresponds approxi¬ 
mately to the value investigated in Fig. 10-I-c. 

C 

When is equal to zero, the poles are located at point 5 and 
indicate the presence of two convergent and two divergent modes. The 

Q 

value or gain of m Increases from zero, at the poles, to infinity along 

^ C 

the direction of the arrows. As m is increased to the value at point 3 

M 

(the break point on the root locus and corresponding to point 3 in Fig. 10B), 
an unstable oscillation starts while the two other roots move along the 

arrows and remain convergent modes. Continuing to increase the gain of 

C 

to larger positive values along point 6, the oscillation increases 
in frequency and instability. The two other poles remain convergences 

Q 

regardless of the positive gain of m . 

C 

The portion of the c.g. travel at which m is negative may be 

* ^ C 

investigated by referring to Fig. 10-IV-c. Again at zero , there are 
the aperiodic modes, two convergences and two divergences (point 5). As 

C 

the gain of is increased negatively, it is seen that a stable oscillation 
starts (point 8) and one divergent mode becomes convergent. However, 
the predominating mode is the ever-increasing divergence (point 7), as 

c c 

m is increased negatively so that the overall mode of motion for ^m 
increasing negatively is an ever-increasing divergence. 

The open-loop response in this case is a serious instability. If 
the roots calculated for the case with the original c.g, location are 
considered (Table I). it is seen that at jx = . 05 and V = 16 mph, there is 
an unstable oscillation with a period of 8.57 seconds and a time to 
double amplitude of only a little over two seconds. In other words, the 
oscillation would double amplitude in about one-quarter of the cycle. 
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4. Long Chord Duct /i = 0.10 V = 34 mph 

Referring to the boundaries of the quartic (Fig. 11A); it is seen 

C 

that m changes approximately 0.0042 per inch of c.g. movement while 

C ^ 

m changes 0,00014 per inch. 
a 

Starting with the center of gravity lowered six inches from the ori¬ 
ginal location, it is noted that the c.g. is on the positive side (point 11) 

of the C = 0 line. As the c.g. is raised toward the original location, the 

C C 

C coefficient changes sign (point 12) and m decreases while m 

fi a 

increases. Continuing to raise the c.g. further, the D coefficient changes 
sign (point 13). When the center of gravity is raised six inches (point 14), 

c 

the line RLBP has been crossed and is negative. If the c.g. is raised a 

few more inches (point 15) the E coefficient changes sign. 

Referring to the stability diagram (Fig. 1 IB);the points again 

are located principally in the region of unstable oscillations. The unstable 

oscillations cease if the c.g. is raised approximately six inches to point 

14 at which location the RLBP line has been crossed and the dynamic 

instability takes the form of a divergence. If the c.g. is raised a few 
C 

more inches, m continues to increase negatively and at point 15, E 
V- 

changes slgn^so that now the configuration is even statically unstable. 

For a closer look at the characteristic modes of motion reference 

is again made to the appropriate root locus plots. From the stability 
C 

diagram note that m ranges in value from approximately 0.001 to 0.004. 

a 

The root locus plots for both these values are shown in Figs. 11-I-d 
and 11-I-e. 

Fig. 11-I-d represents the locus when moving vertically up the 

c c 

line m equal to 0.001 \n the stability diagram (Fig. 1 IB) when m 
« H 

varies from zero to positive infinity. The initial position of the poles, 

C C 

when m equals 0.001 and m equals zero, indicates that two con- 
a p, 

vergences (point 16) and an unstable oscillation (point. 17) exist. As 
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Q 

the gain of m is increased, (equivalent to moving upward along 
C ^ 

m equal to 0.001 in the stability diagram) by moving in the direction 
a 

of the arrows in the root locus diagram, two poles always remain con¬ 
vergences, one with ever increasing value (point 18) and the other 
with decreasing value and a limit at the zero (point 19). The other two 
poles (point 17) continue as unstable oscillations with ever-increasing 

Q 

frequency aqd instability (point 20) as m is increased. 

C ^ 

For m equal to 0.004 (Fig. 11-I-e), the picture is much the same, 
® C 

except that for m equal to zero, the four poles represent two conver- 
P- Q 

gences (point 21) and two divergences (point 22). As the gain of m^ 

is increased the convergences behave as before. However, the two poles 

representing the divergences (point 22) move toward each other and become 

equal roots at point 23, the root locus break point. This point (point 23) 

Q 

is also shown on the stability diagram. As the gain of m^ is increased 
further, the two roots split as shown (point 24) and become unstable 
oscillations of increasing frequency and instability. 

In order to cover the remaining center of gravity locations, (point 

c 

14 and point 15 in Fig. 11B), the root locus plots for m equal to the same 

C » ” 

values, but varying from zero to negative infinity, are analyzed 

(Fig. 11-IV-d and Fig. 11-IV-e). 

Q 

In Fig. 11-IV-d, the four poles at equal to zero represent 
two convergences (point 25) and an unstable oscillation (point 26). As 

c 

the gain of is increased negatively, the two convergent modes 
(point 25) move toward each other in the left half plane to the break point 

Q 

(point 27). As m is increased more negatively, they break away from the 
real axis and become a complex pair representing a stable oscillation 
of increasing frequency and decreasing damping ratio. In the meantime, 

c 

as m is being increased from zero to negative infinity, the unstable 
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oscillation (point 26) is reducing frequency to the point where the complex 
pair of roots (point 26) become divergences, or equal real roots by joining 
on the real axis, (point 28). One divergent pole formerly at point 28 
moves across the origin and becomes a convergence with a limiting value 
at the zero (point 29) in the left half plane. The other divergent pole, 
formerly at point 28, moves to the right, continuing as a divergence of 
ever-increa3_ing instability (point 30). 

Q 

For m equal to 0.004 (Fig. 11-IV-e), the poles behave in much the 
a 

same fashion with the exception that the unstable oscillation at point 26 
becomes two unequal, positive, real roots, representing divergences when 

c 

m equals zero. These poles move in much the same way as the posi- 

M Q 

tive real roots (point 28) move as m^ increases negatively. Note that in 
both cases there is always an instability present, either in the form of 

c 

a divergence or an unstable oscillation for the range of m from zero 
to negative infinity. 

Again, the open-loop response for this case is a serious insta¬ 
bility. Considering the roots calculated for the case of the original c.g. 
location ( Table III), for /i = 0.10 and V = 34 mph, two convergences and 
an unstable oscillation exist. The predominating mode of motion for this 
case is the unstable oscillation with a ten second period and a two second 
time to double amplitude. Therefore, the aircraft doubles amplitude in 
approximately 1/5 of a cycle. 


5,. Long Chord Duct Artificial Stabilization A^ded jlx = 0.10 V= 34mph 
Arbitrarily choosing to investigate the effects of artificial stabi¬ 
lization on the dynamic stability associated with the vertical center of 
gravity located at the original location, an instability exists in the form 
of an unstable oscillation (Fig. 12, points 140 and 150). Taking the case 
of pure attitude feedback (Fig. 12a); one pole cancels with a zero at the 
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origin (T equals infinity), the remaining poles represent two convergences 
(point 141 and point 142) and the unstable oscillation (point 140). As K 
is increased, the unstable oscillation (point 140) becomes neutral (point 
146) and finally stable, with the oscillation reaching a limiting value at 
the zeros (point 143). The two poles at point 141 and point 142 join at 
point 144 and break away from the negative real axis to become a stable 
oscillation of increasing frequency (point 145) as K is increased. 

For the case of pure rate feedback (Fig. 12b), a.zero is located at 
the origin (T equals zero), and the other poles and zeros are located in 
the same position as in Fig. 12a. The unstable oscillation (point 150) 
and the two convergent poles (points 151 and 152) still exist as before. 

However, as the gain K is increased, the poles of the unstable 
oscillation (point 150) move into the left half plane (point 156) and be¬ 
come a stable oscillation with a limiting value at the zeros (point 154). 
Notice that they move in a more direct line contrasted to the former case 
where the loci loop around (point 146) when moving toward the zero (143). 

To investigate intermediate values of T (pure rate feedback with 
various degrees of lag),, refer to Fig. 12a. Locating a zero at the origin 
w<th a pole very close to it but in the left half plane, the loci are 
essentially unchanged. As this pole moves to the left or away from origin 
(T decreasing), the loci between point 140 and point 143 tend to describe 
paths that become more and more direct until they appear exactly like 
the loci between point 150 and 154 when the moving pole reaches negative 
infinity. The other loci for point 141 and 142 remain similar except that 
the break point (point 144) continues to move to the left as the moving 
pole passes the position point 142. The pole at point 142, in this case, 
always moves to the right and terminates at the zero located at the origin. 

From this analysis, it is evident that a pure rate feedback with 
various degrees of lag can be used to stabilize the system and exert some 
control on the frequency and damping ratio of this particular configuration. 
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6. Short Chord Duct ju = 0.05 V = 18 mph 

c 

Referring to the boundaries of the quartic (Fig. 13A); m 
changes approximately 0.002 per inch of vertical c.g. movement while 

Q 

m changes very slightly, say 0.00001 per inch movement. All center 
® C 

of gravity locations have positive m and are located to the left (positive 

a 

side) of the C = 0 line so that the C coefficient in the stability quartic 
is alv/ays positive for this case. However, if the c.g. is raised almost 
six inches above the original location (point 34), the E coefficient changes 
sign so that a dynamic instability appears in the form of a divergence, 
due to the static instability. Raising the c.g. slightly further (point 

Q 

35), m^ becomes negative and an additional increment of upward travel 
causes the D coefficient to change sign as well (point 35). 

Referring to the stability diagram (Fig. 13B); the center of gravity 
located six inches below the original location lies in a region of unstable 
oscillations (point 31). As the center of gravity is raised past the original 

Q 

c.g. location (point 32) ( m decreasing) it crosses the line R = 0 (point 
33), thereby ending the unstable oscillation,and passes through a small 
portion of the stable region before crossing the E = 0 line (point 34) 

Q 

into the divergent region. If the c.g. is raised further, changes 
sign (point 35). 

O 

To detefmine the characteristic modes of motion for this case, 
reference is again made to the appropriate root locus plot. From the stability 

Q 

diagram, m is in the range of approximately 0.0002. The closest 

a Q 

value in the root locus plots is m equal to 0.0005. Analyzing the 

a 

plot (Fig. 13-I-b), which closely resembles the character of the 

Q 

motion for the exact value of m , the following statements can be made: 

C ^ C 

When m equals 0.0005 and m equals zero, one pole (point 36) 

O' ii 

indicates a convergence, the second (point 37) a divergence and the 


' } 


mmn 


remaining two (point 38) a stable oscillation. This is approximately 
the situation existing at point 35 on the stability diagram (Fig. 13B) 

Q 

for the c.g. location six inches above the original. As m is increased 

M 

from zero to positive infinity, the pole at point 36 continues as an increas¬ 
ing convergence (point 48), i.e. , this transient mode tends to die out 
faster. The divergent pole at point 37 moves from the right half plane, 
across the origin (point 34), becoming a convergence in the left half plane 

Q 

and reaching a limiting value at the zero (point 49) when m is increased 
to positive infinity. 

With respect to the two poles at point 38; they increase in fre¬ 
quency and move from the stable left half plane across the imaginary axis 
of neutral oscillations (point 33) into the region of unstable oscillations 
in the right half plane of the root locus diagram. 

It is interesting to compare the appropriate root locus plot to 
the equivalent locations on the stability diagram (Fig. 13B). Notice in 
the root locus diagram (Fig. 13-I-b) the divergent pole at point 37 becomes 
a zero root at point 34 which means that the E coefficient must equal zero 
at that point. 

On the stability diagram.point 34 is located by the intersection 

Q 

of the value of = 0.0002 and the E = 0 line. By closely inspecting 

the stability diagram, it is evident that further increases in ^m result 

\l 

in moving upward across the E = 0 line (point 34) into a region of stable 
oscillations, between point 34 and point 33, on the stability diagram. 

c 

Additional increase of results in moving across the R = 0 
line (point 33) into the region of unstable oscillations. 

On the root locus diagram this means that the pole at point 37 
becomes a convergence (point 34) before the stable oscillation (point 38) 

becomes unstable at point 33. Further increases in on the root 

M 

locus diagram cause the unstable oscillation to become more severe. 
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To cover the remaining c.g. locations which are raised sufficiently 

c 

to make m negative (Fig 13B), the root locus diagram Fig._13-IV-b 

M Q Q 

may be investigated. The value of m is still 0.0005, but m varies 

ot 

from zero to negative infinity on this diagram. The poles located for 

Q 

m equal to zero represent a convergence (point 39) , a divergence (point 
M C 

45) and a stable oscillation (point 4 0 and point 41). As the gain of m^ 

is increased negatively, the poles at point 4 0 and point 41 start decreasing 

frequency as they approach the negative real axis along the direction of 

the arrows Upon reaching the entry.point (point 47) the complex pair 

Q 

becomes tv o equal, real, negative roots. As m continues to be increased 

p 

negatively, the two equal roots at point 47 split and move along the negative 
real axis in opposite directions with one terminating at the zero (point 43) 

Q 

when is equal to negative infinity. Meanwhile, the other pole moves 
from 47 to the left on the real axis where it joins with the original con¬ 
vergent pole (formerly at point 39) at the root locus break point (point 42}. 

C * 

At point 42 they become real ( equal,negative roots and as continues 
to increase negatively they break away from the real axis at point 42 
to become a complex pair, indicating an oscillation of increasing 
frequency (point 44). 

In the meantime, there is always a divergent pole (point 45) 

Q 

which becomes increasingly unstable (point 46) as m^is increased 
negatively. 

Again, reference is made to Table II for the roots of the quartic for 
the short chord duct at/i=0.05,V = 18 mph and the c.g. located in the 
original position. The inodes of motion are two convergences and an 
unstable oscillation with approximately a ten second period and a seven 
second time to double amplitude. This open-loop response,though unstable, 

c 

seems to be more reasonable than the others. In fact, if m is decreased 
by raising the c.g. , the configuration is located in the small stable region 
shown on the stability diagram between point 33 and point 34 in Fig. 13B. 
Actually, this stable region is so small that it is probably insignificant. 



7. Short Chord Duct 


V = 36 mph 


/i = 0.10 

Referri.no to the boundaries of the stability quartic (Fig. 14Aj; 

c 

m changes approximately 0.003 per inch of vertical c.g. movement 
M Q 

while rn changes 0.0001 per inch of movement. All the center of 
a 

gravity locations appear in the divergent region where the coefficients C, 

D and E are negative, except for the lowest location (point 54) where only 
C and E are negative. The highest c.g. location has such a large negative 

Q 

value of m that it is not on the figure. 

M 

Starting in the region labeled'Unstable oscillationd 1 and decreasing 

Q 

m^ along the line of centers of gravity, the E coefficient changes sign 

at point 52. Slightly below point 52 the axis is crossed and ^m becomes 

C ^ 

negative. By increasing m negatively the C coefficient changes 

M Q Q 

sign at point 53 and the location of m and m .for the c.g. located six 

/i a 1 

inches below the originals finally reached (point 54). Raising the c.g. 

Q 

still further, the D coefficient changes sign at point 53 and continues 
to increase negativelv as the higher c.g. locations are approached 
(point 56). 

Referring to the stability diagram (Fig. 14B)i all the c.g. locations 
fall in the region of dynamic instability. The points are below,or on the 
negative side of,the E = 0 line. 

The characteristic modes of motion for these center of gravity 
locations can be obtained from the root locus diagram. 

c c 

Varying m^ negatively at m^ equal to 0.001 approximately 
covers the region in which the centers of gravity are located. In Fig. 

Q 

14-IV-c when equals zero, one pole is located slightly to the right 
of the origin and represents a divergence (point 60). Another pole is located 
in the left half plane and represents a convergence (point 61). The other 
two poles (point 62) are a complex pair located in the right half plane and 
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represent an unstable oscillation. As the gain of is increased 
negatively, the pole at point 61 moves to the right, and has a limiting 

Q 

value at the zero (point 63) when is infinite. This pole (point 61) 
remains in the stable region or left half plane and always represents a 

Q 

convergence which takes place with decreasing rapidity as m is 

• M 

increased negatively 

The two poles at point 62 move along the direction of the arrows 

q 

(point 65) as m is increased negatively. The lightly damped unstable 
M 

oscillation (point 62) increases frequency and becomes a neutral oscillation 
(point 64) and then an unstable oscillation with ever-increasing frequency 
(point 65). 

The predominating mode of motion for this quadrant is the divergent 
pole (point 60) located in the right half plane, or unstable region, when 

q 

m equals zero. This pole moves to the right (point 66) and causes 
M 0 

the transient motion of the machine to become more unstable as m 

is increased negatively. 

If the center of gravity is lowered considerably more than six inches, 
0 

through some design change, m becomes positive. Investigating the 

M 0 0 

root locus using the extrapolated value of m » 0.001 and m varying 

a /i 

from zero to positive infinity, the root locus appears as Fig. 14-I-c. 

C 

The poles for m^equal to zero, indicate a convergent mode at point 67, 
an unstable oscillation at point 69, and the fourth pole (point 68) close 
to the origin is actually in the right half plane (shown by close analysis 

Q 

of Fig. 14-I-h or Fig. 14C-b for m increased positively to 0.001 
C a 

while m equals zero) and represents a divergence. As the gain of 

C ^ 

m is increased positively from zero, the pole representing a divergence 
(point 68) moves to the left toward the origin and into the left half plane 
to become a convergence (point 71). The pole at point 67 always moves 
to the left (point 70) and remains in the left half plane as a convergence. 
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The complex pair (point 69) remains unstable and represents an oscil¬ 
latory mode of increasing frequency and instabilityjpoint 72). 

C C 

For the case of the actual c.g. locations, m positive and m 

a /x 

negative, as shown on the stability diagram, the open loop response is 
a serious instability. 

C C 

Taking the specific case of the values of m and m at the 
C * 

original c.g. location (where m is admittedly very unstable), the 
roots of the characteristic equation (Table II) for /x = 0.10 V = 36 mph 
clearly define the modes of motion. The convergence and stable oscilr 
lation are strongly masked by an instability in the form of a divergence 
that doubles in amplitude in less than a second ( T - 0.8068 seconds). 


8. Short Chord Duct Artificial Stabilization Added /x = 0.10 V = 36 

For the short chord duct at /x = 0.10, there is a dynamic instability 
in the form of a divergence. 

Analyzing the limit cases of pure rate feedback (T = 0) and pure 
attitude feedback ( T = oo), the root loci appear as in Fig. 15. The pole 
at point 120 represents the instability due to the divergence and the 
zeros are located for increasing K from zero to infinity. When T is equal 
to infinity, the pole and the zero located at the origin cancel each other. 
When T is equal to zero, the pole is located at negative infinity and the 
zero remains fixed at the origin. 

In the case of pure attitude feedback (Fig. 15a), where the pole 
cancels the zero at the origin (T equal to infinity) the remaining poles 
represent a stable oscillation ( point 115), a convergence (point 117), 
and the divergence (point 120). As K is increased the stable oscillation 
(point 115) continues, but increases in frequency (point 116). The poles 
at point 117 and point 120 move toward each other until they become 
equal roots at the root locus break point (point 119) in the right half 
plane. As K is increased further, the poles break away from the positive 
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real axis ( point 119) and become an unstable oscillation, a neutral 
oscillation (point 121), and finally a stable oscillation with a iim.ting 
value at the zeros (point 118) in the left half plane. 

For the case of pure rate feedback (Fig. 15b), in which a zero 
(point 136) is located at the origin (T equal to zero), the same modes of 
motion prevail as before when K is equal to zero. However, in this case, 
as K is increased, the poles at point 130 move toward the break point 
(point 131) and split so that one moves toward point 132 and the other 
toward point 133. The pole reaching point 133 joins the pole formerly at 
point 135 to break away from the real axis into the stable oscillation 
which reaches a limiting value at the zeros (point 134). Meanwhile, the 
pole (point 137) representing the divergence can only become a zero root 
at point 136 when K is increased to infinity. 

In order, to study the effects of various degrees of lag, an extra 
pole is situated in the left half piano, very clo.^e to the zero (point 136). 
As K is increased, it will join with the pole at point 135 and break away as 
a stable oscillation terminating at the zeros (point 134). The poles at 
point 130 remain as stable oscillations which move away from the real 
axis and continue as stable oscillations of ever-increasing frequency. 

As this extra pole is moved to the left, away from the zero (point 
136), the former description of the loci remains much the same except as 
it moves past the vicinity of point 131 toward point 132. When the 
extra pole moves in this direction it pulls the ioci of point 130 closer and 
closer to the real axis until, at some large negative location, the poles 
(point 130) join on the real axis similarly to that shown at point 131. 

The pole located at point 137 always moves to the zero (point 136) 
and the basic path is not influenced by the extra pole moving in the left 
half plane. 
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From this analysis it is possibl i to note that pure attitude feedback 
can stabilize this system (point 121 and point 118) for some value of gain, 
K, whereas pure rate feedback can only yield a zero root (point 136). 

The gyro-bar with either aerodynamic or mechanical damping 
represented by the transfer function for a lagged rate feedback, 


JCd 


1 + Td ' 

should not be referred to as a method of providing rate and attitude feedback. 
These devices do not sense attitude and therefore provide only rate feedback 
with a phase lag. While it may be true that a lagged rate feedback provides a 
frequency and phase characteristic similar to that of combined rate and 
attitude feedback, this is only true within a certain frequency range so 
that, strictly speaking, there really is no direct equivalence. 

For this reason the gyro-bar devices may be able to diminish the 
degree of dynamic instability of an aircraft in forward flight associated with 
a divergence due to a static instability, (E < 0), but they can never completely 
stabilize the system. This is exactly the case for the short chord duct at 
ji equal to 0.10 and the e g. in the original location. Even at infinite gain, 

K, the divergent pole (point 137) can only be made a zero root (point 136). 

An understanding of the cnaracteristic modes of motion for configura- 

Q 

tions of the Flying Platform with other combinations of values for m and 

C ^ 

m^ not discussed is possible by reference to the other root loci provided 

for all the quadrants of each stability diagram. 
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Generally, the unstablli z«d, seven-foot diameter, ducted rotor 
Flying Platform exhibits a degree of control-fixed dynamic instability 
with either the long or short chord duct for all vertical center of. gravity 
locations. 

In hovering, both duct configurations are dynamically unstable 

c * 

due to the large positive values of m . This instability is characterized 
by an unstable oscillation for all center of gravity locations considered. 

In forward flight, both duct configurations exhibit dynamic 
instability characterized by unstable oscillations for the lower c.g. 
locations and divergences for the higher c.g. locations. 

For any of the configurations, artificial stabilization, using 
a rate feedback with a phase lag as provided by a gyro-oar, is able to 
correct a dynamic instability when it Is of the form of an unstable 
oscillation. When the instab iity is characterized by a divergence the 
lagged rate feedback device cannot completely overcome the instability. 
However, a combination of rate and attitude or pure attitude feedback 
can stabilize this divergence. 

The commonly used gyro-bars with aerodynamic or mechanical 
damping can only provide ratp fpedhack with a phase lag and never provide 
any of the attitude feedback required to completely correct the instability 
associated with a positive real root of a characteristic equation. However, 
lagged rate feedback may certainly help to minimize this type of an 
instability even though it can never completely overcome it. 

The large tilt angles oredicted for the faster trimmed flight speeds 
will become a problem for eithe r duct configuration. For the ducted rotors 
studied, these excessive tilt angles cannot be attained because the pilot 
has insufficient control moment input to trim the large moments generated 
by the aircraft for the higher speeds. 

In general, for given c.g. locations the moments produced by the 
short chord duct were less than the long chord duct, but the moment could 
be reduced for eithr. duct by raising the center of gravity. 
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Any attempt to increase the attainable forward speed by compromising 
between these duct designs or center of gravity locations does not 
produce significant or desirable results. Reducing the moment by raising the 

c.g. causes the divergent mode to appear (point where the duct becomes 

aCM 

statically unstable or^— = E = 0) at successively lower flight speeds. 


y 
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TABLE I 


CHARACTERISTIC EQUATIONS OF THE LONG C fORD DUCT 7 1 

Center of gravity located, in the original position (0.42c) 

T - Time to double or half amplitude P - Period of oscillation 


Stability Cubic ^overing [l = V = 0 t = 0.3864 

A. + 0.31272 A + 0.046425 A + 0.37872 * 0 


Roots 

A = - 0.81973 

T = O .32666 seconds 

Roots 

a 2 = 0.2535 + 0.630691 
= 0.2535 - 0.630691 

T = I.O 563 seconds. 

P = 3*8495 seconds 


Stability Quartic /i = 0.05 V = l 6 mph T = 0.3944 

A 4 + 0 .53207 A 3 - 0.00007544 A 2 + 0.055276 A + 0.00868 = 0 

Roots 

A = - 0.13976 

T = 1.956 seconds 

Roots 

A 2 = - 0.63364 

T = 0.43135 seconds 

Roots 

* 3 = 0.12066 + 0.288871 
x 4 = 0.12066 - 0.288871 

T = 2.265 seconds 

P = 8.5786 seconds 


Stability Quartic fJ. = 0.10 V = 34 mph T = 0.3809 

A 4 + 0.62572 x 3 - 0.0031 A 2 + 0.0282 x + 0.010023 = 0 

Roots 

* 4 = - 0.21140 

T - 1.2486 seconds 

Roots 

^2 = - O .66036 

T = 0.39973 seconds 

1 Roots , 

*3 = 0.12302 + 0.238051 
a 4 = 0.12302 - 0.238051 

T = 2.1457 seconds 1 

P =10.054 seconds 















TABLE II 
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CHARACTERISTIC EQUATIONS OF THE SHORT CHORD DUCT 
Center of gravity located in the original position (0.42c) 


T - Time to double or half amplitude P - Period of oscillation 


Stability Cubic Hovering ^ = V = 0 r = 0.33926 

x 3 + o. 34493 x 2 + 0 . 014260 x + 0.073105 = 0 

Roots 

X = -0-55575 

T = 0.42304 seconds 

Roots 

X 2 = 0.1054 + 0.346841 
* X 3 = 0.1054 - 0 .346841 

T = 2.2306 seconds 

P = 6.1459 seconds 


Stability Quartic ^ = 0.05 V = l8 mph T = O.35I 

x 4 + 0.45072 x 3 + o.o4o44 x 2 + 0.02577 x + 0.0013625 * 0 

Roots 

x l = -0.05511 

T = 4.4l4 seconds - 

Roots 

X 2 -0.4684? 

T = 0.5192 seconds 

Roots 

+ 0.03643 + 0.226821 
+ 0.03643 - 0.226821 

T = 6.677 seconds 

P = 9-7231 seconds 


' Stability Quartic M= 0.10 V = 36 mph T = 0.359 

x + 0.54622x 3 + 0.00029x z ' - 0.058851A - 0.006944 « 0 

Foots 

Ay 4* 0.30838 I 

T = C.8068 seconds 

Roots 

a 2 = - 0.13518 

T = 1.8484 seconds 

Roots 

X 3= - ®. 35971 + 0.192851 

X 4= - O.35971 - 0.192851 

T * O.69163 seconds 

P = 11.6964 seconds 









TABLE III 

OVERALL VALUE Of LIABILITY DERIVATIVES 
FOP. VARIOUS C. 0. LOCATIONS 


FEET DOWN 


LONG CHORD DUCT 


N 


ORIGINAL 


FEET UP 


o 

II 

0 

0.1772 

0 

0. Io 04 

0 

0.1435 

0 

0.1266 

■ 

O.OOPl 

0.0510 

0.002? 

0.0324 

0.0023 

O.OI36 

0.0024 

- 0.0051 

m 

0.0005 

0.0506 

0.0014 

0.0250 

0.0022 

- 0.0005 

0.0032 

- 0.0266 


SHORT CHORD DUC r 


C.G. LOCATION C.G. LOCATION C.G. LOCATION C.G. LOCATION 


FEET’ DOWN ORIGINAL I 0.5 FEET UP 



=0 

0 

0.0465 

0 

'o.o 4 o? 

0 

0.0339 

0 

0.0276 


0.00005 

0.0239 

0.0001 

0 0121 

0.0002 

0.00013 

0.00025 

- 0.0119 

A 

.10 

0.0011 

-0.0099 

0.0017 

-0.0272 

0.0025 

- 0.0446 

0.0032 

- 0.0621 








































































FIGURE 1 FLYING PLATFORM 







:ZE DUCT PROFILES FOR WIND TUNNEL MODEL TESTS 







All notation shown is positive except flight path 
angle to horizon , 7 . 


FIGURE 3 BODY AXIS SYSTEM 











FIGURE 4 BODY AXIS INERTIA TERMS - CURVED FLIGHT PATH 








Velocity Ratio 



FIGURE 5 TYPICAL TIME RESPONSE TO A SMALL INPUT IN HOVERING 










Velocity Ratio 



FIGURE 6 TYPICAL TIME RESPONSE TO A SMALL INPUT IN FORWARD FLIGHT 




















FIGURE 7 ROOT LOCUS S.C.D. AND L.C.D. HOVER 
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c.g. locations 
0 0.5 ft. down 

O original 

0 0.5 ft. up 


1.0 f*. up 


LONG CHORD DUCT : 

STABLE OSCILLATIONS 

STARt AT C = 0.0003. 
m 

UNSTABJJE OSCILLATIONS 
START AT C^ = 0.0056. 

_ M 


SHORT CHORD DUCT: 

STABLE OSCILLATIONS 

START AT C = 0.000075 
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M 

UNSTABLE OSCILLATIONS 
START AT C m = 0.0024 


FIGURE 8 STABILITY - DIAGRAM FOR HOVER C > 0. 
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FIGURE 9 ROOT LOCUS S.C.D. HOVER-ARTIFICIAL STABILIZATION ADDED 
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STABILITY DIAGRAM 


V = 16 MPH 


FIGURE IOB 
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FIGURE 10 C ROOT LOCUS - L.C.D. - ji » 0.05 V = 16 NFH 



































C.D 




































FIGURE 10 - IV 


ROOT LOCUS 


L „C„ D. /i = 0.05 


V = 16 MPH 







0.5 ft. down 


original 


BOUNDARIES OF THE STABILITY i 


FIGURE 11A 
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UNSTABLE 
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O 0.5 FEET DOWN 
O ORIGINAL 
O 0.5 FEET UP 
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DIVERGENCE 


UNSTABLE REGION 


STABILITY DIAGRAM 


L.C.D. = 0.10 V = 34 MPH 

FIGURE IIB 
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BOUNDARIES OF THE STABILITY QUARTIC S.G. D. fi = 0.05 V = I8MPH 


FIGURE 13A 
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STABILITY DIAGRAM 


FIGURE I3B 
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BOUNDARIES OF THE STABILITY QUARTIC S.C.D. a = 0.10 V = 36 


FIGURE 14A 
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STABILITY DIAGRAM S.C.D. \l «0.I0 V * 36 MPH 

FIGURE 14 B 
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FIGURE 14 - II 
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E VALUATION OF STABILITY DERIVATIVES 

The stability derivatives were evaluated from data obtained by 
wind tunnel tests conducted at the David Taylor Model Basin. A two foot 
diameter,. 2/7 scale, model of a ducted-rotor Flying Platform was tested 
in the DTMB low-speed wind tunnel (Ref. 2). From the large variety of 
test runs made, portions of the data incorporating different duct designs 
were selected for the purpose of obtaining stability derivatives. The duct 
designs chosen were: one of relatively long chord length with a large 
leading edge radius, the other of shorter chord length and smaller leading 
edge radius (Fig. 2). The latter was chosen because it represented a con¬ 
siderable design change from the large leading edge, long chord duct. 

The tests were conducted in the large settling chamber of the 
wind tunnel. The balance was designed so that forces were measured 
perpendicular and parallel to the rotor shaft and the moments were 
measured about the supporting structure. 

» 

The moment was arbitrarily transferred to the so-called quarter 
chord point oi the duct. The forces and moments so measured were T, 

H and M c/4 (Fig. 3). 

All runs were made at constant velocity ratios as the angle of attack 
of the duct was varied. The tests were conducted principally at only 
three different velocity ratios, making accurate determination of derivatives 
with respect to p, somewhat difficult. 

When summing all the forces and moments acting on the machine, it 
was necessary to include the drag force and moment caused by the pilot 
in addition to the forces and moments measured in the wind tunnel, (Fig, 3) 
This arag force, directed opposite and parallel to the relative wind, was 



i »• 





taken as acting through an assumed drag center of the pilot. In this report, 


the drag center for the typical pilot is located 3- feet above the pilot's 


Various methods are available to determine the drag coefficient or 
equivalent flat plate area of the pilot. A report published on the subject of 
airloads on human beings (Ref. 9), determines the drag coefficient of 
var'o; s men in standing positions. For the specific case of a six-foot , 

200 pound, clothed man, standing erect, with arms and feet together, the 
equivalent flat plate area was determined as approximately nine and one 
half square feet. No attempt was made to modify this value due to the 
pilot stance with arms extended and feet apart, or to compensate for ad¬ 
ditional components, such as the cage around the pilot, instrumentation, 
fuel tanks, pilot helmet and equipment. It was felt that attempts to make 
accurate modifications on the drag of the pilot, in any way, would not 
be worthwhile unless investigations were made of the air flow to which 
the pilot is exposed. Therefore, for the purposes of this report, the 
value of the equivalent flat plate area of the pilot is assumed to be 
nine and one half square feet. 

The original center of gravity position for the complete machine 
at normal gross weight was arbitrarily located at the pilot's stand-on 
platform or 0.42c above the top of the duct. Additional vertical center 
of gravity positions were selected above and below this original reference 
point (Fig. 16). All moments were resolved about these center of gravity 
locations and the stability derivatives evaluated for each case. (Table III). 

The important requirement in the evaluation of the stability deri¬ 
vatives is that they be calculated at each initial, steady-state, equili¬ 
brium trim condition of the aircraft. For this steady flight condition, the 
sum of all the aerodynamic forces in the drag direction must be equal to 
zero and the sum of the forces in the lift direction must equal to the gross 
weight of the aircraft. 






4 














Fig. 20 Equilibrium Diagram 


In Fig. 2 0 where or is a small angle, the summation of forces 
in the drag direction is: 


This equation may be solved to determine the variation of the 
equilibrium angles of attack with forward speed. 

Actually, this was accomplished graphically by plotting the 
horizontal force versus angle of attack for the various forward speeds 
(Fig. 21). 











Values of equilibrium angles plotted versus speed are shown for 
both ducts in Fig. 18 on page 117. 

For vertical equilibrium in the initial steady-state level flight 
condition, (when the sum of the forces in the lift direction equals the 
gross weight) it is necessary to determine the variation of rotor angular 
velocity with forward speed from: 

L = W= (Vc r a*■*«.<)/> 7rR 2 (nR) Z < 64 ) 

where R * 3.5 feet, W* 550 pounds and C TT , and a are obtained 

H T 

from the wind tunnel data for the particular equilibrium angle and velocity 
ratio. 

Then, for vertical equilibrium the actual flight speed may be 
determined from: 

y= a QfR 

' (65) 

Actual values of fiR at the various velocity or tip speed ratios 
are listed for both ducts in Table IV. 
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TABLE IV 

ROTOR TIP SPEED (it/sec) versus yi 


Velocity Ratio 

j o 

ii 

fi= 0.05 

M ~ 0.10 

*2R Short Chord Duct 

555 

537 

525 

| ftR Long Chord Duct 

--- L 

433 

478 

495 


When evaluating the stability derivatives, the forces and moments 
were put into coefficient form and plotted versus the angle of attack 
(measured in radians) and ji (dimensionless velocity ratio). 

Typical sketches in the evaluation of the derivatives are shown 
in the following Figures. 



The short vertical lines, on the constant p curves, indicate 
the previously calculated equilibrium angles for each duct at each 






velocity ratio. The stability derivative must be evaluated at the 
equilibrium points by taking the slope of the curves at these points. 

Q 

Actual values of T versus a are listed in Table V. 
a 


TABLE V 


T versus u. 
a 


Velocity Ratio 


T Short Chord Duct 
a 


T Long Chord Duct 0 


ll = 0.05 

fi = 0,*1 

0.003 

0.011 

0.0047 

0.0133 


To obtain T^ it is.necessary to cross plot the former Figure 


of C„ veisus a. 
T 

For C T : 


-— 

Q 

- = o 

_a * -0.1 

a - -0 .7 

a = -0. 




3 




a = -0. 

/L = equ 

_ 

librium a 

t - - 

1 

\ a = - 

L,- 


F f ig. 23 Sketch of C^ versus p, 


* • * 











The stability derivative T^ is obtained by taking the slope at 
the equilibrium angle for the particular pL. 

Occasionally it was necessary to extrapolate in order to obtain 

the hovering values {p. * 0) of derivatives taken with respect to p. 

This can be avoided if more values at zero alpha are obtained at the 

higher p or by additional data at lower velocity ratios. 

C 

Values of T obtained from the data are listed in Table VI* 


TABLE VI 

c: 

T versus u 


Velocity Ratio 

0 

M = 0.05 

p = 0,10 

C 

T^ Short Chord Duct 

0 

0.044 

0. 016 

c 

T Long Chord Duct 

0 

-0.004 

-0.044 


For C H : 


O' 



Equilibrium angles of attack for /x = 0.15 were usually so 
large that they were not shown on the figures. 




Values of ~H are listed in Table VII. 
a 


TABLE VII 


C versus /i 


Velocity Ratio 


i* = 0 /i = 0.05 \ il = 0.10 


0^ Short Chord Duct q 

a 

Long Chord Duct q 

L_r._^_ 


0.0005 | 0.0015 


Cress plotting for C. 


a - 0, -0.1, -0.2 ^ o- = 



A = equilibrium a 


0.05 0.10 0.15 

V- 

Fig. 25 Sketch of C„ versus n 
H 


Values of C u are listed in Table VIII. 

h m 




DPI 

NS 









TABLE VIII 


C versus p 
H 


Velocity Ratio 


H Short Chord Duct I 0.044 


n=0 

/l = 0.05 1 

/!= 0. 

10 

0.044 

0.078 

0.108 

1 

CO 

r-H 

r—t 

o 

I 0.118 

0.131 

5 


For C 


A.= equilibrium a 




- 1.2 - 0.6 
o- (radians) 


Fig. 26 Sketch of C versus a 
m 


Values for C are listed in Table DC. 
m 

a 


TABLE DC 


C versus a 
m 

O’ 


Velocity Ratio 


m Short Chord Duct I 0 
a 


m Long Chord Duct I n 
a ! 






‘V 


*. .*• . ’I 


.. Vv 


Ml*. jk. ... 


.• f ’ 














I 



Cross plotting f. 


A = equilibrium o- 


Fig. 27 Sketch of C 


Values for m versus |iare listed in Table X 


TABLE X 
m versus u 

It r 


Velocity Ratio 


m Short Chord Duct 


0.055(ext.) 0.03 


m Long Chord Duct 


The evaluation of the derivative C. 


may be obtained by 


3[d(A0). 
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Values for H 


ii, z 


H/R 


d (A9) 
, and r are listed 


using the appropriate values 
in Table XI. 


of 



TABLE XI 


C versus \x 

H d U8) 


Velocity Ratio 

1! i 

o 

= 0.05 

(i = 0.10 

C H ,, ..Short Chord 
d(A8) Duct 

0.000275 

0.000486 

0.000673 

c 

H,, . Long Chord 

d ( A 0) _ 

v 7 Duct 

0.000735 

0.000735 

0.000848 


In the equations of motion the drag terms (Terms 5, 6 and 7 
on page 19) are developed as a result of the drag of the pilot and may be 
evaluated directly. Their form in the equations may also be seen in 
Fig. 20 where the coefficient form of this force (acting parallel to the 
relative wind) is // C^' /2 (Eq. 13). 

Using the small angle assumptions, the pilot drag force, 

2 

u C„ l /2 is resolved into components along the body axes. The derivatives 
D 

are taken with respect to | 4 .and a and must be evaluated at the initial 

condition, a - u. and a = a . 

2 ° 

X component: il C n ' 

2 D 
2 ^ , 

Z component: \i a 


Taking the derivatives: 


X axis: C D ' 
*2 



2 


_a_ 

8m 

_a 

da 


(jl ) = llC D ' 

(M 2 ) = 0 


o 


= “o C D 









J 



< Tlm% 




'Aia; a (. ■ 


2 15L 


-■ "fl _3 W * u 2 C-' = M o 


The forces in coefficient form are obtained by multiplying these 
terms by their respective perturbations* 


H-force contribution: ^ C • A/* (Term 5 page 19) 


Thrust force contributions: C D * A/x (Term 6 page 19) 


2_ Cp' Atr (Term 7 page 19) 




-J»r 
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The remaining stability derivative to be evaluated is the over-all 

c 

pitch damping derivative m^ (Eq. .21g), 


c = c 

d(^6) d(A 9) 


-in 

R 


C » 

d(AQ) 


(21g) 


The total damping moment in pitch for the unstabilized Flying 
Platform includes contributions from three sources, namely: 

1. The pitch damping moment due to the change in H-force with 


pitch rate, hL.C 


H 


. 3.(±e) 


R XA d ( a 9) 

2, The pitch damping moment of the rotor due to a pitch rate. 

3. The pitch damping moment of the duct due to a pitch ate. 

Contributions (2) and (3) taken together are equai to the first term 
C 

m .,. in Eq. 21g. 


d(A<?) 


Contribution (1) is obtained by solving Eq. 68 for il 


d(A9) 


C 


H d(*e) 


~ r . i±L , / 

/a R Hy 


The values of H 


zpj, D , R and t are all known. 


Then 


the pitch damping moment contribution due to change of H-force with 
pitch rate is: 


( 68 ) 


ikt . • d(nO) 

R d(*6) 


(Term j. 7 page 21) 


Contribution (2) is obtained by using the expression developed 
in the Appendix of Ref. 4 for the moment generated by the rotor due to 
a pitch rate. 




- -ji tin) = - (&e) dl 

dt 

~(^e)(p7rffv e )(jS. 2 i-c*) 


C. 

dfA0> 

*OTO t « 


- _ a ^ c -Pd R 

/6 W 


where is the average chord of the rotor blade. 

Contribution (3); the pitch damping moment of the duct due to a 
pitch rate may be derived in various ways. 

The moment produced by this contribution is equal to the time 
rate of change of angular momentum. For a constant pitch rate, this is 
equal to the time rate of change of moment of inertia of the slipstream 
at the duct exit which is equal to the rate of mass flow times the square 
of the radius of gyration for the duct exit area about the pitching axis. 
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c = 

*oro* 


This.gives: 


(Afl) 


- pitch rate about 
Y axis at top of 
duct. 

- duct exit area 

- duct chord length 

- slipstream velocity 
at duct exit. 











umm 


where; 


pnR~ V^, = rate of mass flow 

o 

R" 2 

4 + c - square of the raaius of gyration for the duct exit area 


The duct exit velocity is equal to (See Appendix Ilf): 






Therefore the pitch damping moment in coefficient form is; 


r. (*9)f*R £ nR\ c? r/?' 

~~ "p»R*(CIR)*R IT* 

r - (<£a0){C? PtR 3 /? f / r*1 

~ ~ : y* —x 

sm /? J 


The pitch damping derivative is; 

where w is equal to the aircraft relative density factor m/prrR 2 R. 

For the long chord duct when R = 3.5 ft. and c =1.75 ft, c 2 /R 2 = 1/4 
The expression for contribution (3) is; 


d(LQ) 


I£lL 

2 ur 
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The sum of the three moment contributions gives the total theoretical 
pitch damping moment used in this Report. 

It is recognized that this theoretical pitch damping moment is probably 
too large, since it is not very likely that all three contributions provide their 
full theoretical value for the cases of the relatively short chord ^ucted rotors 
studied. In any event, this total value appears to be small compared to the 
amount required tc produce desirable stability characteristics. Therefore, 
the stability characteristics together with the effects produced by additional 
damping were investigated and are liscussed in this Report. 
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APPENDIX III 


The expression for duct exit velocity may be developed from 

the so-called simple momentum theory of ducted-propellers by utilizing 

the common assumption that the velocity of the air at the duct exit 

is equal to the final slipstream velocity (V . = V ). 

exit oo 

The analysis Is restricted to the inviscid, incompressible, 
static case and utilizes the ideal ducted propeller combination with an 
‘•actuator c isk" and a simple cylindrical shroud. 


The velocity V is the same on both sides of the actuator disk 
because of the requirement of continuity. 

The total thriist is: 


Utilizing the assumption of no slipstream contraction for a 
ducted pi speller Ap = A^, * A . 

By continuity: 
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